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ABSTRACT 

Gamma-ray bursts (GRBs) are emitted by relativistic ejecta from powerful cosmic explosions. Their light 
curves suggest that the 7-ray emission occurs at early stages of the ejecta expansion, well before it decelerates 
in the ambient medium. If so, the launched 7-ray front must overtake the ejecta and sweep the ambient medium 
outward. As a result a gap is opened between the ejecta and the medium that surfs the radiation front ahead. 
Effectively, the ejecta moves in a cavity until it reaches a radius -R gap ~ lO^-Eg^cm where E is the isotropic 
energy of the GRB. At R = i? gap the gap is closed, a blast wave forms and collects the medium swept by radiation. 
Further development of the blast wave is strongly affected by the leading radiation front: the front plays the role 
of a precursor where the medium is loaded with e ± pairs and preaccelerated just ahead of the blast. It impacts the 
emission from the blast at R < R\ oa a = 5i? gap (the early afterglow). A spectacular observational effect results: 
GRB afterglows should start in optical/UV and evolve fast ( < rnin ) to a normal X-ray afterglow. The early optical 
emission observed in GRB 990123 may be explained in this way. The impact of the front is especially strong if the 
ambient medium is a wind from a massive progenitor of the GRB. In this case three phenomena are predicted: (1) 
The ejecta decelerates at R < i?i oa d producing a lot of soft radiation. (2) The light curve of soft emission peaks at 
t P cak ~ 40(1 + -z)i?54 2 (r o j/100)~ 2 s where r o j is the Lorentz factor of the ejecta. Given measured redshift z and 
£pc a k, one finds r o j. (3) The GRB acquires a spectral break at 5 — 50 MeV because harder photons are absorbed 
by radiation scattered in the wind. A measurement of the break position will determine the wind density. 

Subject headings: Cosmology: miscellaneous — gamma-rays: bursts — radiation mechanisms: nonthermal — 
scattering — shock waves 
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1. INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) are explosions of 
huge energy ~ 10 52 — 10 54 ergs (see Piran 1999 for a review). 
The relativistic ejecta of the explosion produces the observed 7- 
ray pulse with duration of a few seconds that propagates ahead 
of the ejecta and interacts with the ambient medium first, before 
the blast wave driven by the ejecta. 

Madau & Thompson (2000) and Thomson & Madau (2000, 
hereafter TM) pointed out that the 7-ray pulse can preacceler- 
ate the ambient medium to a high Lorentz factor. Even more 
importantly, the pulse-medium interaction is accompanied by 
runaway loading of e pairs (TM). The interaction occurs in- 
side the thin radiation front where the primary photons scatter 
off the medium and turn into e ± via 7 — 7 reaction; the cre- 
ated e ± increase the medium opacity, do more scattering, and 
next generations of e ± are loaded in a runaway manner. TM 
also pointed out that the e ± loading and preacceleration ahead 
of the blast wave should modify the GRB afterglows (see also 
recent paper by Meszaios, Ramirez-Ruiz, & Rees 2001). 

In the present paper we develop an accurate model for the ra- 
diation front and assess its impact on the blast wave. The prop- 
agating front is self-similar and its non-linear structure (the rise 
of density and velocity across the front) admits a simple de- 
scription. We find the medium parameters behind the front and 
identify the range of radii where the impact on the medium is 
strong and hence the ensuing blast wave is strongly affected. 
Especially interesting effects are found for GRBs with massive 
progenitors, leading to spectacular observational phenomena. 

In this paper we assume that the radiation front is formed 
early inside the ejecta (the so-called "internal" scenario which 



easily explains the fast variability of GRBs, see Piran 1999). 
The front energy remains constant with radius i.e. we neglect 
additional radiation from the ensuing blast wave (afterglow) 
compared to the prompt radiation from the ejecta. In a separate 
paper (Beloborodov A.M., in preparation) we study the oppo- 
site case where the ejecta emission is negligible and the 7-ray 
front is created by the blast wave itself (the "external" scenario). 

In §§ 2 and 3 a detailed formulation of the problem and ba- 
sic equations are given. Numerical solution is presented in § 4. 
In § 5 we develop an analytical model that explains the front 
structure and reproduces the numerical results with good accu- 
racy. The backreaction of the GRB-medium interaction on the 
prompt 7-rays is studied in § 6. The sweeping of the medium 
by radiation and the front evolution with radius are studied in 
§ 7. In § 8 we compute the blast wave dynamics in the preac- 
celerated environment and evaluate its emission. 

2. FORMULATION OF THE PROBLEM 

2. 1 . Basic parameters of the front 

GRB produces a thin shell ("front") of collimated radiation 
with bolometric flux F(m) and spectrum F e (xu) where e = 
hv/m e c 2 and vd is the Lagrangian coordinate in the moving 
shell, < n7 < A. Here A is the front thickness and A/c is the 
observed duration of the burst. The front propagates through 
the ambient medium with velocity c. The medium interac- 
tion with the front is convenient to view in the nj-coordinate: 
medium "enters" the A-shell at w = 0, passes through it and 
goes out at vj = A with new density and velocity. Radiation 
scattered by the medium is decollimated and also streams to- 
ward large za, being absorbed by the primary beam. 
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RADIATION FRONT OF GAMMA-RAY BURSTS 



The scattering of GRB radiation can have a strong impact 
on the medium if each electron scatters many photons during 
its passage through the A-shell. The photons "kicked out" by 
the electron from the collimated 7-ray beam can be converted 
into e , so a large number of scatterings would imply a large 
number of pairs created per one ambient electron. The main 
contribution to pair production comes from photons with e ~ 1 
(see § 5), and their density is n p h ~ F/m e c 3 . The electron 
scatters many photons if its "free path" A = 1 /n p hCT (the dif- 
ference 8m between successive scatterings) is smaller than the 
front width, 2 

Tfh C 

X = — ^— = 4.64 x lO 6 ^!^ 1 cm < A. (1) 

The radiation flux is F = L/AttR 2 where R is the distance 
from the center of the explosion and L is the isotropic lumi- 
nosity of the GRB. The total energy of the radiation pulse is 
E = (A/c)L and the condition ([!]) can be rewritten as 

R<R X = (^y^ 8.0 XlO^f cm. (2) 

Beside A there is another important length-scale in the front 
— the typical 5zu the scattered photons pass before they get 
absorbed by the primary radiation. This "photon free path" 
(hereafter denoted A 77 ) far exceeds A. It implies that pair cre- 
ation occurs at larger w i.e. substantially lags behind scat- 
tering. As we show in § 5, the runway pair loading starts at 
m = a w ^/AA 77 w 30A. Efficient pair loading in the front 
requires A > 30A which implies a tighter constraint on radius: 
R < i?ioad = R\/VS0. At radii larger than i?i oa( j there is 
neither pair loading nor acceleration of the medium. 

-Rioad should be compared with the deceleration radius of the 
GRB ejecta. The standard R dcc ~ 10 16 - 10 17 cm > i? load 
holds if the ejecta decelerates in a normal interstellar medium 
(ISM). However, if the GRB has a massive progenitor, its am- 
bient medium is the wind from the progenitor; then i?dcc ~ 
10 14 — 10 16 cm depends on the mass loss M and the velocity w 
of the wind, e.g. R dcc ~ 10 15 cm for M = 10" 5 M Q yr" 1 and 
w = 10 s cm s" 1 expected for a Wolf-Rayet progenitor (Cheva- 
lier & Li 1999). The impact of the radiation front is especially 
strong in the latter scenario where i?dcc < -Rioad- In particular, 
as shown in § 8, the very value of i?dec is changed. 

The front model we construct makes use of the fact that 
the scattering optical depth of the ambient medium tr = 
<TTRpo/m p is very small. The GRB pulse-medium interaction 
occurs in a specific regime: each electron of the medium expe- 
riences a lot of scattering while GRB photons have a low proba- 
bility of scattering. (Pair loading increases the optical depth; in 
the calculations we assume that the medium remains transpar- 
ent and discuss the conditions under which it becomes opaque 
in § 6.1.) Another important feature of GRBs is that the radi- 
ation density strongly dominates over the rest-mass of the am- 
bient medium, F/c ^> poc 2 . The primary radiation dominates 
over the scattered radiation, e , and magnetic field in the front. 

When the backreaction on the radiation pulse is negligible on 
time-scales < R/c, the propagating front is quasi-steady. It can 
be formalized as follows. Let us define 

w = ct — R, (3) 

2 Standard notation is used throughout the paper: a magnitude Y measured in 



where t is the time passed since the beginning of the explosion. 
Then we have R = ct and w = at the leading boundary 
of the front, and to = A at the trailing boundary. Now let us 
change variables (t,R) — ► (t : m). That the front is quasi-steady 
means that the medium parameters are functions of vj <C R and 
t w R/c is a slowly changing parameter. The front gradually 
changes when its radius R increases. We aim to construct a 
model for the front structure, i.e. determine the medium den- 
sity and velocity as functions of w. We will show that the front 
evolves in a self-similar manner: with increasing R, its struc- 
ture is given by same unique functions of dimensionless coor- 
dinate £ = w/X. 

The front is thin (A <C R) and its quasi-steady structure 
(formed on time-scales <C R/c) can be described in plane- 
parallel geometry. The collimation of GRB radiation is very 
strong, A8 < 0.01, and to a first approximation it is perfectly 
collimated. Asymmetry of the explosion does not affect the 
calculations unless the ejecta are beamed within an angle less 
than AO. Possible inhomogeneity of the ambient medium is not 
considered in this paper [in contrast, Dermer & Bottcher (2000) 
discussed the impact of the 7-ray s on circumstellar clouds]. 



2.2. Particle collectivization and the cold approximation 

It is reasonable to assume that the loaded e ± share imme- 
diately their net momentum with the medium. Even though 
the Coulomb collisions are extremely inefficient, there are two 
other mechanisms of momentum exchange: 

(1) The created e ± form a stream interacting with the 
medium via beam instability. The instability time-scale is of or- 
der cjpj 1 where o; p i = (47rn e e 2 /m e ) 1 / 2 is the plasma frequency 
and n e is the electron density. 

(2) In the presence of transverse magnetic field B the pairs 
gyrate around the field lines frozen into the medium on the Lar- 
mor time, uj^ 1 = rn e c/Be. The net momentum of is thus 
communicated to the medium. 

The first mechanism should always work because its time- 
scale is shorter than the Compton cooling time of e , tc = 
(3m e c 2 /8cttF), and also much shorter than the time of 
medium dynamics across the front. The magnetic coupling 
may dominate if lob > ^pi which requires B 2 /An > n e m e c 2 . 
When the medium accelerates one should substitute the rest- 
frame magnetic field, density, and flux in these estimates. 
The acceleration results in compression (Madau & Thompson 
2000), the density and magnetic field are amplified, and the 
coupling becomes even stronger. 

The coupling passes the net momentum of injected to the 
medium and maintains their distribution approximately isotro- 
pic in the medium rest frame. It does not ensure however that 
the also share their energy with other particles. One can 
distinguish between two possible situations: (1) Partial col- 
lectivization: e ± injected with a Lorentz factor j e (measured 
in the medium rest frame) quickly get isotropized but preserve 
7 e ; their subsequent cooling is controlled by Compton losses. 
This is the case with magnetic coupling. (2) Complete collec- 
tivization: all particles share their energy instantaneously and 
maintain a Maxwellian distribution in the medium rest frame. 
It might be the case if collective modes provide sufficient cou- 
pling. 

We will show that the majority of e^- are loaded with mod- 
erately relativistic energies, cool efficiently (even with partial 

units and divided by 10 k is denoted as Y^. 
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collectivization), and remain at subrelativistic energies. Par- 
ticles created at w are Compton cooled much faster than the 
medium moves to zu + a where next generation of hot pairs is 
created. Therefore the bulk of pair-creating radiation at zu + a 
has been scattered by cooled particles. To the first approxima- 
tion, the medium can be considered as a cold plasma with a 
bulk velocity f3 found from momentum conservation. We here- 
after use this "cold" approximation since it greatly simplifies 
the calculations; its validity is checked in § 4. 

3. BASIC EQUATIONS 

In this section we give the equations of a steady radiation 
front in the plane -parallel geometry (see § 2.1). 

3.1. Scattering and pair creation 

Let p, be the cosine of the scattering angle. A primary colli- 
mated photon scattered through p starts to move backward with 
respect to the A-shell with velocity dzu/dt = c(l — fi) and its 
ru-coordinate grows. The scattering at < zu' < zu determines 
the intensity of scattered radiation at zu, 



IsC (Ml ^Sc) 



1-p 



2tt dp e 



(4) 



Here dzv'/(l — p) = cdt is the length element along the scat- 
tered ray, da /dp is Compton cross-section (see Appendix), n 
is the electron/positron density, and j3 is the medium velocity in 
units of c. F e , n, and f3 are taken at the location of scattering, 
zu' . The photon energies before and after the scattering, e and 
e sc , are related by 



e(l-/3) 



1- /3 M +(l- M ) e / 7 ' 



(5) 



P 2 ) 



-1/2 



is the Lorentz factor of the scattering 



where 7 = (1 
medium. 

The scattered radiation that propagates from zu' to zu is at- 
tenuated by 7 — 7 absorption. This is accounted for by the 
exponential factor in equation (Q) where t 77 is the 7 — 7 optical 
depth, 



vdn7. 



(6) 



The opacity dominated by the primary collimated ra- 

diation (the scattered radiation has much smaller density, see 
§2.1). A scattered photon (e sc ,p) can interact with primary 
photons e that are above the threshold 



Cthr 



(1 - mK 



(7) 



The cross-section for interaction with e Jj e t hr is c 77 ~ O.Iitt. 
The r 77 can be viewed as the product of (1 — /i)er 77 and 
the column density of primary photons above the threshold, 
~ s(F tthr /m e c 3 ) where s — c(t — t') = (zu — zu')/(l — p) is 
the path passed by the scattered photon. 

The exact expression for «; 77 (/i, e sc ) is given in equation (Q) 
of Appendix. In numerical examples we will consider a homo- 
geneous primary radiation pulse i.e. assume that the spectrum 
F e does not depend on 137. Then the 7 — 7 opacity is homoge- 
neous across the A-shell and r 77 (n7, zu') = (zu — zu')k^. 



The pair creation rate at given zu is the rate of 7 — 7 interac- 
tion between I sc {p>, £sc) and the primary beam F e , 



n + (zu) = 2n J dp J de 

rGJ f 1 p 

J dm' J dfi J de 



^^(1-^)^-^.(8) 



Here we made use of equation (Q). 

3.2. Continuity equation 

Let rii and n e be the density of background ions and elec- 
trons, and let 2n+ be the density of created pairs. The total 
electron density of the medium n = n e + 2n + and its veloc- 
ity v — (3c satisfy the continuity equation. For a plane-parallel 
front the continuity equation reads 



dn d(nv) 
~dt + dR 



2n + — 2h a 



(9) 



where h + and n aim are the local rates of pair creation and anni- 
hilation, respectively. The annihilation rate r'i ann = (3/8) (1 — 
/3 2 )n^_(TTC is many orders of magnitude smaller than h+ and 
hereafter we neglect annihilation. 

Since both n and v are functions of zu = ct — R only, we 
have d/dt = cd/dzu and d/dR = —d/dzu, and rewrite the 
continuity equation as 



(10) 



The immediate consequence of this equation is that the magni- 
tude 



n(l-(3) 



(11) 



would conserve in the absence of pair creation and hence the 
compression of accelerated medium is (1 — /3) _1 (see also 
Madau & Thompson 2000). In particular, for the background 
electrons and ions we have 

n* e ee n e (l - (3) = n , n* = m(l - (3) = n i0 . (12) 

Here no and ni Q are the electron and ion densities prior to the 
interaction with the front. 

The mass density of the medium is (we neglect the additional 
mass associated with the plasma internal energy: the cold ap- 
proximation) 



Po 
1-/3 



1 



n [i e m v 



Me 



n Q m p ' 
(13) 



We neglected the small contribution (~ m e /m p ) of the back- 
ground electrons to pq. The is the ion mass and p, e is the 
medium mass (in units of m p ) per electron: p e = 1 for hydro- 
gen and [i e — 2 for helium or heavier ions. The ratio n* /no 
shows the number of e^- loaded per one background electron. 

The cross-section for Compton scattering is inversely propor- 
tional to the squared mass of the scatter, so only are efficient 
scatters. The average mass per one scatter is 



P 
n 



(14) 



The initial m* 
loading. 



p e m p can decrease to m e as a result of pair 
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3.3. Momentum conservation 

The law of momentum conservation reads (neglecting the 
pressure forces: the cold approximation) 



d(vjp) d(v 2 jp) 



dt 



dR 



— P± + -P SC ) 



(15) 



where P± is the momentum deposited by pair creation per unit 
volume per unit time and P sc is the momentum deposited by 
photon scattering off the medium. We rewrite this equation as 



(16) 



The scattering passes momentum from the beamed radiation to 
the medium with rate 



7 4 \ n* f F e 
1 — J — / de— / dcr(e - /j,e s 



(17) 



The factor 1 — 7 4 /7sat accounts for a finite collimation angle of 
the primary radiation (see eq. ^ of Appendix). Assuming that 
the radiation is emitted by the ejecta with Lorentz factor T e j at 
R = i? om , we have 7 sat = T cj (R/R cnl ) at a radius R. 
The momentum deposited by pair creation is given by 



P± 



dm' / d/i / de s 

I! J -I 



22 n * e r„ (18) 

em e c z q/x 



Here p± (/x, e sc ) is the average momentum of the e ± pair created 
when a scattered photon (fi, e sc ) gets absorbed, 



P± 



(19) 



The numerical factor x ~ 1 is given in equation (g) of Appen- 
dix. 

3.4. Thermal balance 

The continuity and momentum equations allow one to com- 
pute the dynamics of the medium in the cold approximation. 
When we know the dynamics of the cold medium, we can eval- 
uate its temperature from the thermal balance; it will allow us to 
check the consistency of the cold approximation. The thermal 
balance in the medium rest frame reads 



d(uV) 



dV 



dt -P^+^ c 



,d (nV) 
dt 



+ {C + -C-)V. (20) 



Here u is internal energy density of the medium (including rest 
mass of e^), p is pressure, r )- lu ^{m)m e c 2 is the mean energy of 
injected e ± , V is volume per barion, and dt = dt/7; all these 
magnitudes are measured in the rest frame of the medium. 

The terms C ± are the rates of Compton heating/cooling. 
Both depend on the particle energy distribution in the medium 
rest frame. Given the uncertainty of this distribution we re- 
place it by (5-function at a mean Lorentz factor 7 e and estimate 
roughly 



C+ -C~ 



{Phl - Pill) nVa T F T 



(21) 



where 7 e = 7c corresponds to Compton equilibrium and Ft is 
the flux of (primary) radiation that scatters in Thomson regime; 



this flux is measured in the medium rest frame and it is approx- 
imately 



F T — Fr(e < ckn 



1-/3 7(1+1) 

7-7-3, £KN ~ • (22) 

1 + P le 



Here £kn is the typical energy above which the scattering oc- 
curs in the Klein-Nishina regime. Fr(e < ekn) is the primary 
flux with e < 6kn> measured in the lab frame. 

The ions carry a small fraction of the thermal energy (even if 
they manage to share the energy with e ± , their density rijCn 
as soon as pair creation begins) and hence u ~ j e m e c 2 n. Note 
that nV = n* and V = 7(1 — (3)Vo where Vq is volume 
per barion prior to the interaction with the front. Substituting 
these relations into equation ( |2C| ) and taking into account that 
dva/dt = 07(1 — (3) we get after simple algebra 



d7 e (p/hm e c 
dm 



' d [7(1 _^ )]+ (^j-7e)dn* 



7(1-/3) dm 1 ' v r n ' n* dm 

, 4 (7c ~ 7e) °"tPt 
3 7(1 + (3) m e c 3 



(23) 



We estimate the pressure interpolating between nonrelativistic 

p/hm e c 2 — (2/3)(7e — 1) and relativistic p/nm e c 2 = 7 e /3 
limits, 



P 



7<: 



1 



nm e r 



37e 



kT 

m e c 2 



e. 



(24) 



Here we introduced a temperature T. For a Maxwellian plasma 
T is related to pressure by p = hkT in both non-relativistic and 
relativistic cases. For a non-Maxwellian distribution, T is an 
effective temperature defined by p = nkT. 

Once we know ^(m) and n*(m) from the dynamic "cold" 
solution we can find Ft and 7i n j(w) (see Appendix). Then we 
can solve numerically equation ( |23| ) and find 7 e (o7) and Q(m). 

4. NUMERICAL RESULTS 

In this section we construct a numerical model of the front. 
Here we assume that the ambient medium is hydrogen (/i e = 1, 
see eq. |l3[); the extension to 1 < \i e < 2 is simple (§5). 

We integrate the ordinary differential equations ( |lo| ) and ( [To] ) 
with the boundary conditions (3 = and n = rio at m — 0. 
At each step dm we know the radiation scattered at previous 
steps (smaller m) and find the local pair creation rate from 
equation (Q) and the rate of momentum injection from equa- 
tions (|l7|) and (|l8|). After getting the dynamic solution n(m) 
and 7(07) we integrate the thermal balance equation (^3]) with 
the boundary condition 7 e (0) = 1 and find je^m). 

The input of the calculations is the GRB spectrum F £ (m) and 
the output is the front structure n(m)/no, j3{m), and ^^m). 
For numerical illustration we take a radiation pulse with con- 
stant spectrum 



Fie- Q S e<l, 
F, = I Fe~ Q2 , K e < e br , 
0, e > e hr . 



(25) 



Such a spectral shape is observed in GRBs with at\ ~ ±0.5 
and a?2 ~ 1.5 ± 0.5 (Preece et al. 2000). In numerical exam- 
ples we fix a\ = and assume a.% > 1. Then the problem has 
a well defined solution in the limit ebr — ► 00 ■ The finiteness of 
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£br ^ 1 causes a break in pair loading at w = m^r (see §§ 5.3 
and 6.2). In the examples below we take = 10 2 . 

The solution is a function of the dimensionless coordinate 3 
£ = zu/X where A is given by equation ([!]) with the total flux 



F 



(a 2 - ai)Fi 
(1 - ai)(a 2 - 1) 



(26) 



4.1. Dynamical structure of the front 



Figure 1 shows the computed front structure in the case of 
«2 = 1.5 and 7 sat = 10 3 (the solution does not depend on 7 S , 
until 7 approaches 7 sat , cf. 



eq. 



17) 



Near the leading bound- 
=! 1 and pair loading pro- 
30. At £acc « 10 2 the 



ary of the front, the medium has 7 » 
ceeds exponentially on scale of £i oa d 
medium accelerates. 

The acceleration length can be understood in simple terms. 
When a portion dm of the radiation pulse overtakes an electron 
with 7 w 1, it passes momentum dp m 0.2(F<tt / c 2 )dvj (here 
0.2 is a Klein-Nishina correction). Hence dp/d£ « 0.2m e c 
and the medium acceleration length is £ acc « 5m* /m e where 



yields an estimate 



(rio/n*) is mass per scatter (see eqs. [l4| and|l3[). This 



5(m p /m e ) exp(-£ acc /30), 



(27) 



i.e. £ acc ~ 10 2 . The estimate neglects the additional accelera- 
tion due to P± which is approximately equal to P sc at £ ~ £ acc 
(see Fig. 2). More exact formulae are derived in § 5. 

The pair loading continues in the accelerated zone £ > £ acc 
and m* further decreases. Therefore the medium accelerates 
very fast, 7 w (£/£ aC c) 3 , until n* /no reaches m p /m e and m* 
saturates at m e ; afterward 7 cx £ 3 / 2 . 

The accelerating medium scatters radiation through smaller 
angles, /j, w /3 — > 1, and pair loading slows down: d 2 n*/d£ 2 
becomes negative at £ w £ aC c- The decrease in pair produc- 
tion is caused by the growth of the threshold for 7 — 7 inter- 
action, e t hr oc (1 — /i)" 1 (eq. 0). The scattering by medium 
with a relativistic 7 produces photons with 1 — fj, as I/27 2 , 
and hence e t hr oc 7 2 oc (c;/Cacc) 6 - The 7 — 7 opacity seen 
by the beamed scattered photons, 



becomes very 

low and they travel almost freely across the front with a free- 



thr 



path 



(C/Cacc) • For instance photons with 



e = 1 scattered by medium with 7 = 2 (at £ = 1.3£ acc ) are 
absorbed only at £ ~ 10 4 , i.e. absorption strongly lags behind 
scattering. It implies that photons scattered at £ <; £ acc con- 
trol the pair loading in the whole relativistic zone of the front. 
The pair loading thus decouples from the medium dynamics at 
£ 3* £ aC c- As explained in § 5, the bulk of radiation scattered at 
£ > £ sc « 1.4(e br /100) 1 / 6 (where 7 > 7sc = v^/4) is not 
absorbed at all and escapes the front. 

At £± sa 10 3 the e 1 * 1 density exceeds the density of the am- 



bient electrons b 
Then equations ( 



the factor m p /m e 
and (O) yield 



and to* saturates at m P 



^ = ~ 1 



Tsat / 



7 dn 
n* d£ 



(28) 



a good approximation at 7 3> 1 (see Fig. 2). In the absence of 
pair loading (dlnn*/dln£ <C 1), 7 would tend to the asymp- 
totics 7 = £/2. However before 7 can reach any asymptotics it 
saturates. The saturation happens at £ = £ c ^ 10 4 . 
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Here we took (3 sa 1 and calculated P sc with Thomson cross- 
section (medium with 7 ^> 1 scatters the GRB radiation in 
Thomson regime). We neglected P± compared to P sc , which is 

3 The ^-coordinate has the meaning of dimensionless fluence of the burst, £ = (ctt /m e c 3 )Fvu. The computed n(£) are 7(g) are also the exact solution for bursts 
with arbitrary light curves F\ (to) (but with a fixed spectral shape) once £ is defined as £ = (or/meC 3 ) Fdza'. 



FIG. 1. — Structure of the radiation front for ot\ = 0, 02 = 1-5, tbr = 10 . 
Top: Dynamic structure. Dashed and .vo/irf curves show i.*(O/ Tl ar >d 7(C) 
where £ = ro/A. So/i'rf vertical line shows the boundary of the ion-free zone, 
5c, whe r e 7 reach es 7sat- Dotted curves show the analytical model of § 5 (see 
eqs. f!9|p5[|52[p3| ). Bottom: Thermal structure. Solid curve shows the mean 
kinetic energy or particles in the medium rest frame. The other two curves 
display 7i n j — 1 and 7c — 1 (cf. the text). 



Our steady dynamic problem becomes inconsistent when 7 
saturates. The assumption that the front has the speed of light 
and the medium passes through it with dm /At — 1 — (3 be- 
comes wrong. Instead the medium gets stuck in the front: it has 
reached the velocity /3 sat such that the net flux of GRB radiation 
vanishes in the medium rest frame. The /3 sat is determined by 
the angular spread of the radiation and represents the effective 
velocity of the radiation pulse. Saturation implies that the am- 
bient (ion) medium is trapped in the pulse and cannot penetrate 
the zone £ > £ c — this zone is ion-free. (More exactly, the ions 
cannot penetrate fj > max{£ c , £ mix } where £ mix ~ 7~ 2 i?/A, 
see § 7.1.) The trapped ions accumulate and surf the pulse. 
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FIG. 2. — Momentum deposition rate in the front shown in Figure 1. Solid 
and dashed curves display P BC and P± (see eqs. |l6||li| ). 

Radiation scattered by the medium in the process of its ac- 
celeration partially propagates to the ion-free zone of the front 
and produces e ± there owing to 7 — 7 reaction with the pri- 
mary radiation. A steady h + (zu) is established throughout the 
whole front on a relatively short time-scale ~ wj (l—/i)c where 
1—fi ?S (l/27g C ) ~ 0.1 represents the typical collimation angle 
of the scattered radiation that produces pairs. The pairs created 
in the ion-free zone acquire the saturated Lorentz factor, stay 
almost static in the m— coordinate, and accumulate. 

The constructed plane-parallel model can be applied to the 
expanding spherical front as long as the time-scales involved 
do not exceed R/c. In particular, the time i acc (7) of medium 
acceleration to a given 7(137) should be smaller than R/c. To 
estimate £ acc assume the most favorable conditions for acceler- 
ation: to* = to 6 and dlnn*/dln£ <C 1; then equation ( j28| ) 
reads d-f/dt = (c/A)(l - /3)/(l + /3) » (c/4 7 2 A) and gives 
£ aC c ~ 7 3 Vc- From i acc < R/c we find 7 < 7 max where 



7d 



1/3 



-1/3, 



'15 



1/3 
'53 • 



(29) 



Both 7 sat and 7 max evolve as the spherical front expands. If 
7max < 7sat the boundary of the ion-free zone £ c is determined 
by 7 = 7 max rather than 7 = 7 sat . We discuss the front evo- 
lution in more detail in § 7 and derive £ c (R) there. Here it is 
worth to emphasize that the front structure shown in Figure 1 at 
£ < £ c does not change with radius, i.e. the front is self-similar. 
Since A oc R 2 the front is just "stretched" in tn-coordinate as 
R grows. Note also that the front trailing boundary £a = A/ A 
in Figure 1 moves to the left with increasing R and becomes 
smaller than £ c at some radius R c found in § 7. At R > R c 
there is no ion-free zone and the whole front is described by the 
self-similar solution. 

4.2. Thermal structure of the front 

In the bottom panel of Figure 1 one can see two peaks of 7 e 
at £ « 70 and £ 4 x 10 3 . They correspond to the beginning 



and the end of the medium acceleration. This unusual temper- 
ature profile can be understood from equation ( j23| ) which we 
rewrite as 



d7 e 
d£ 



e 



7inj - 7e 



£a 



load 



7e 



(30) 



Here the effective temperature is related to the average 
Lorentz factor j e via equation (]24j), £ acc = (— dln[7(l + 
/3)]/d£) _1 is the acceleration length, £ load = (dlnn*/d£) _1 
is the pair loading length, and 



£c 



3 F 7(l + l)7e 

4 F T 7 e 2 -7c 



7e > 1, 



3 F 7(l+/3) 

4 F T 7 e • 

3 F 7(l+/3) n „ 1 

4 F T 01-131 ' < Je ^ X > 



(31) 



is the length of Compton cooling. 

The initial temperature of the medium is low and it gradu- 
ally rises at small £ owing to injection of pairs with 7; n j ~ 10. 
Already at £ 3 the temperature exceeds the Compton equi- 
librium value 0c ■ Thereafter > ©c and Compton scatter- 
ing cools the medium rather than heats. At £ £ioad ~ 30 
the pair density exceeds that of background electrons and be- 
gins to exponentiate. One could then expect a high heating 
rate, however Compton cooling is very efficient and keeps the 
temperature below m e c 2 . The length of Compton cooling is 
£c ~ (F/Ft)ekn ~ 1. It is much shorter than £i oa d and there- 
fore the cooling competes successfully with the heating. This 
competition is described by equation 



d7e _ 7inj _ 7e 

d£ £i oa d £c 



(32) 



Here we neglected the first (adiabatic) term on the right-hand 
side of equation (|o|) since it is much smaller than the other 
two terms. Equation ( |32| ) shows that 7 e — 1 saturates at ~ 
£c7inj/£ioad ^ 1. This is the first maximum of the tempera- 
ture profile. 

At £ « £ acc w 10 2 the medium begins to accelerate and 
then the relative velocity between the injected e ± stream and 
the medium decreases. Correspondingly, 7j n j, the heating rate, 
and the medium temperature fall down. 

When the medium Lorentz factor reaches 7 ~ 10, the rela- 
tive velocity between the injected e stream and the medium 
vanishes and changes sign. Here 7i n j reaches a minimum. Af- 
terward the e loading tries to decelerate the medium (see also 
§ 5.4). The acceleration by scattering, however, dominates and 
the medium continues to accelerate. Now 7i„j rises again (the 
relative velocity between the injected e ± and the medium again 
increases) and the heating rate and the temperature grow. 

The cold approximation is especially good near the mini- 
mum of 7 e at about 2£ acc . The temperature is also quite low 
at £ acc < £ < 2£ acc where the main scattering occurs (that con- 
trols pair loading in the whole accelerated zone of the front). 

The energy distribution of e around the average 7 e depends 
on details of their thermalization. In the case of partial col- 
lectivization (see § 2.2) the distribution has a tail extending 
from 7 e to 7; n j whose slope is controlled by Compton cool- 
ing. The length-scale for Compton cooling of injected e ± is 
£c(7inj) ~ {F/F T )e K N ~ 1 (it does not depend on 7 inj or 
7 because F^/F ~ £kn in the case of a± = 0). Hence the 
density of pairs with 7 ~ 7; n j is 



dn* 



aload 



(33) 
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i.e. the number density of the high-energy particles is about 
30 times smaller as compared to cooled particles. Also their 
energy density is smaller than that of the cooled component. 

We conclude that the cold approximation is reasonably good. 
Note however that we focus on relatively soft spectra oti > 1 
(in contrast, TM took «2 = 1 as a basic case). The case of hard 
spectra «2 ^ 1 is more complicated because the maximum of 
7e at £ £ acc becomes essentially relativistic and then numeri- 
cal simulations relaxing the cold approximation will be needed. 
A relativistically hot plasma scatters preferentially backward 
(smaller /i); note also that P± then strongly dominates over P sc 
at £ ~ ■face- We expect however that the front structure will not 
change qualitatively for hard spectra, though the values of £i oa d 
and £ acc may change. 

5. ANALYTICAL MODEL 

The medium dynamics in the radiation front can be under- 
stood with a simplified model that we formulate below. In par- 
ticular, we derive the characteristic lengths £i oac i and £ acc , get an 
analytical solution for the front in the non-relativistic (0 < 0.5) 
zone, and evaluate the pair loading rate in the accelerated zone. 

5.1. Formulation 
Let us replace the scattering cross-section by 



= a T S(p - 0)H(e KN - e), 



(34) 



where S is the Dirac function and H is the Heaviside step func- 
tion. Here we have made two approximations: 

1. Assume that radiation scatters with Thomson cross- 
section if e < 6kn and does not scatter at all if e > 6kn, 
where £kn Ss 1 is the energy above which the Klein- 
Nishina corrections reduce the scattering and subsequent 
pair creation. We derive in Appendix the effective 



£kn ps 0.47(1 + 0) 
for calculations of n + and P±, and 

e£S«0.7 7 (l + /?) 
for calculations of R r . 



(35) 



(36) 



2. Replace the broad distribution of the scattering angles by 
its average, /i = 0, i.e. assume that the collimated ra- 
diation scatters through 90° (p, = 0) in the medium rest 
frame. Then we also have 



(37) 



1 + 



The scattered photons can interact with primary photons of 
energy e > ethr where the threshold is given by equation (fy. 
In our simplified model equation (Q) reads 



Ethr 



2(l+ff) 

6(1 - 0) ■ 



(38) 



We have ethr > 1 for any e < £kn, the scattered ra- 
diation interacts with the high-energy part of the spectrum, 



F e = Fie ° 2 , e > 1. The 7 — 7 opacity of the power-law 
radiation seen by the scattered photon is (eq. HQ in Appendix) 



K 11 = —■ \ H ( £ thr 



— ^(ethr-ebr), A i"^- 



(39) 



The numerical factor (p(a) can be approximated with high ac- 
curacy as (Svensson 1987) 



(a)«l 2 -<*(l + a)- 5 / 3 . 



(40) 



Hereafter we use notation (j> = (j)(ct2) =0.045 and 0.023 for 
oti =1.5 and 2 respectively. Equation ( ^9| ) is exact for a 
power-law spectrum and inaccuracies appear only when ethr 
approaches the spectral break et> r - Given the opacity, we also 
know the free path of the scattered photons A 77 = kzI, 



. _ Ai /£ t hr\ Q2 
A77 ~ J \~) • 



(41) 



Finally, let us replace the exponential attenuation of the scat- 
tered radiation in equation (Q) by the step function H(l — t 77 ). 
Then equations (§), (111), and (p/Th read 



xe a ^ 1 H{e max -e)H{e-e mia ), (42) 



X( J Q 1 + 0' \1 + 0'J Jo 

xe a -- 1 P± iJ(e max - e)H{e - e min ), (43) 



Psc = 1 - 



iLJ (l+/?)Ai Jo 



fede. 



(44) 



Here f e = F/F 1 = e~ Ql if e < 1 and f e = e"" 2 if e > 1, 
£min is found from the condition e t hr < £br, 



(45) 



and e max is found from the condition r 77 = [w — ru')K~ n < 1, 



M = 



Ai 



(w — ZD 1 ) 



l/a 2 



1+0' 
1-0' 



(46) 



When e max > et,r one should replace the upper limit by £bi - 
This refinement is however not important since eb r is anyway 
far from the scattered peak e ~ 1 that dominates pair loading 
as shown below. 

The formula for p± (eq. [FT] of Appendix) in our simplified 
model reads 



P± 



cB' 



(1 + 0') 0(a 2 -l) 



m e c 1 + 0' e(l-0') 



(47) 



"2 



The e-integrals in equations ((42|), (^3j), and (^) depend on 
the relative positions of e m i n , e max , ckn, and unity. We now 
consider two different zones of the front starting from small w. 
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5.2. Non-relativistic zone ((3 <C 1) 
In the non-relativistic zone we have e m ; n <C £kn 



and equations ( |Tq| ) and ( f42| ) give 

dn 2n + 20 
~~ ~ _ A? 



< 1 < 



dec? 



J / dro'n(tu') / e Q2 -7 e de 

1 -/O ie mm 

n(m') dm' . 



20e" 2 " Ql 



KN 



(48) 



\\(a 2 - a\) 

Here we neglected e m i n compared to ekn- The exact solution of 
equation (^) is the sum of growing and decaying exponentials, 

1/2 



vj /a 



m/a) 



a 



( a 2 



Oil 



\20e" 



(49) 



no ( 

Substituting ol\ = and 6kn = 0.4 we get £i oa d = a/A w 24 
and 33 for a 2 = 1.5 and 2 respectively (here we used Ai/A = 
F/Pi given by eq. ^). The analytical solution is in perfect 
agreement with the numerical results, see Figures 1 and 3. 

The loading length admits an easy interpretation. As seen 
from equation (^8j), scattered photons with e ~ £kn make the 
dominant contribution to 
free-path of these photons, 



Equations Q) and (J38|) give the 



A 77 « 



— C*2 

KN ■ 



(50) 



One can see that a ^/A 77 Ai ss ^/A 77 A. Note that A 77 /A ~ 
200 and 500 for a2 = 1.5 and 2 respectively, i.e. the scattered 
radiation is weakly absorbed in the non-relativistic zone (this is 
a consequence of £kn < Cmax)- When an ambient electron has 
passed a distance m through the front, it has scattered ~ m/X 
photons and a fraction ~ m of these photons have been 
absorbed. Hence one pair is injected per one ambient electron 
when (m/X) x (ro/A 77 ) = 1 which gives the above formula 
for the loading length a. 

We now evaluate the medium acceleration at /3 <C 1. First 



let us calculate the momentum loaded b 
p± from equation into equation ( 



pairs. Substituting 



we get 



P± 



4>m e 



A? 



dm n(m') 



dei 



$(ct2 - 1) 

<j>(a 2 ) 



(a 2 - l)m e 



-KN 



-Oi\—\ 



2\{{a 2 -ai-1) 



■ noa I e 



vj / a 



3 — w/a 



(51) 



This is a perfect approximation if ebr — > oo. 



in ( |51[ ) peaks at the upper limit as e^N " 



The e-integral 
One can take 

e * ~ £ kn/2 as a typical e of scattered photons, then one gets 
the typical energy of absorbed primary photons e a bs ~ X £ thr ~ 
(1 + a 2 1 ) 5/3 (2/e^) ps 20, assumed to be well below e br . Note 
that a 2 — cti — 1 when a 2 — > 1 and a?i = 0, i.e. for hard 
spectra the e-integral in j5l| ) does not have a pronounced peak 
and low-energy photons e -C £kn contribute a lot to P±. Such 
scattered photons interact with very energetic primary photons 
Cabs ~ IOckn/c and then the finiteness of ebr is important: 
£abs > £br is excluded which leads to a reduction of P±. E.g. 
in the case of a 2 = 1.5 and e br = 10 2 the actual P± is sup- 
pressed by a factor of 2 compared to equation (j5l|). 

Equation ( p| ) gives a perfect approximation to the momen- 
tum deposited by scattering. Where /? <§C 1 it yields 



™ e c 2 (egS) 



acc A a i + 1 



Ai ai + 1 2 



L>*>/a 



-vj J a 



(52) 



The medium accelerates according to the momentum equa- 
tion (p"oj). With 7 w 1 and p ~ po this equation reads 



d/3 
dnj 



P+ + P« 



(53) 



Substituting (plh and ( 52 ) and integrating for {3 we get 



2p e m p Xi 



(a 2 - l)e 



KX 



l)Ai V 6 



^—•uj/a 



(a 2 — ai — 1)A 

/ accNQx+l 
| \ fc KNi / e T37/a 

ai + 1 V 



7 /a 



(54) 



Here we used po/no = p e m p (see eq. 13). The non-relativistic 
zone ends when (3 reaches ~ 0.5. Equating j3 = 0.5 and ne- 
glecting the decaying exponential we get the acceleration length 
(with ai = 0, e K N = 0.4, and efgj = 0.7), 



In. 



(p e m p /m e )(\i/a) 
7^?M«/Ai) + 



5 + ln/i e . (55) 



0.7 



Hence £ acc w 5£i oa d at /i e = 1, in full agreement with the 
numerical simulations (Fig. 1 and 3). As one can see from 
equation (p5J), with n e — 2 the result changes only slightly, 
£ acc w 5.7^i oa d. Note that £i oa d does not depend on p e at all. 
Hence the front structure is not sensitive to the chemical com- 
position of the ambient medium. 

5.3. Relativistic zone ( (3 — > 1) 

At m > m acc the medium continues to accelerate relativisti- 
cally. Then 6kn grows (eq. [35] ) and exceeds unity. The integral 
over m' in equations (^) and (^3]) is now taken over two re- 
gions: < m' < mi where e max (m') < 1 and mi < m' < m 
where e max (m') > 1. The boundary w\ is defined by condition 

^max 1^ 



m — mi 



Ai A + A 



l-/3i 



(56) 



This is an implicit equation for vj\ where j3i — (3(mi). One 
can show that m acc < mi <C m when 7(117) ^> 1. From equa- 
tion (E3) we then find 



0c r , n*(m') fl-(3' 



Q(m'), 



(57) 



O(w') 



oci—ol\ 



m < nji, 

lnmin{e max ,e K N}, m' > m x . 



The integral peaks at m' ~ n7i (where e m i n <§C e ma x ~ 1 
and we therefore set e m i n ~ in the expression for Q). De- 
note the integrand as S and evaluate the integral as J dm' S w 
£zj7i>S(ti7i). We have 5 oc n* at ttj' < n7i and a steep decline 
5 oc n*y 2a2 at to' > m x , hence C « (dln?i*/dlnn7 + 
From the numerical results we see that C 1/3. Then we get 



h+ (m) 



C,cn\mi 



(a 2 - ai)Xi(m — mi) 



(58) 
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This formula gives a reasonable approximation to n + at w > 
ro acc (see Fig. 3). 

The approximation e m i n <C 1 used in the derivation of equa- 
tion (|8) breaks when e m i n (ci7 1 ) approaches unity i.e. e t h r 
for scattered photons with e = 1 approaches £br- We define 
a characteristic zu sc such that e m i n (ci7i = ra7 sc ) = 1/2 (i.e. 
ethr = £br/2). The velocity of the scattering medium at w sc is 
given by (see eq. p3h 



7s. 



/Cbr 



if £ b r > 100 



(59) 



At w\ > w sc the scattered peak e ~ 1 is not absorbed at any 
w. The corresponding cut off in pair loading appears at 



nj br 



Al /£br\ Q2 



Cbr\ a2 



/e b r\ Q2 _ /Cbr\ 



(60) 



We conclude that (1) pair loading at any vo is sensitive to the 
medium dynamics at w < w sc only and (2) the extension of 
the pair loading zone is limited by finite £br- 

The simple qualitative picture of pair loading in the relativis- 
tic zone is as follows. The scattering of photons with e ~ 1 
makes dominant contribution to h + at any w < W\, T . Pho- 
tons scattered at a given vo\ > m acc get absorbed at w = 



A 77 where 



i + A 



ft > 0.5. 



(61) 



The scattering in a narrow interval £ acc < £ < 3£ acc controls 
pair loading in the whole relativistic zone £ acc < £ < 10 8 . 
Unfortunately, we do not have any simple analytical solution at 
£acc < £ < 3£ acc . The empirical formulae 



( £ fee ) ' £acc < £ < 3£ a( 

3V3(^) 3/2 , £>3£ acc , 

(f^) 2 , Cacc <£<3£ a 

3(|£), e>3e acc , 



(62) 



(63) 



fit well the numerical results for both a 2 = 1.5 and a 2 = 2 
(see Fig. 1 and 3). Here n* « 0.5/i e e 5 rin ~ 74/j, e no is known 
from § 5.2 (eqs. |49|andp5]). 

5.4. Heating by pair loading 

We now give estimates for the mean energy and momen- 
tum of the injected pairs, first in the lab frame (e±, p±) and 
then in the rest frame of the medium (7i n j, Pinj)- The estimates 
highlight the role of e ± in the heating and acceleration of the 
medium. 

In the lab frame, the energy and momentum of a created e 
pair is dominated by the absorbed primary (collimated) photon, 



m e c 2 



P± 

m e c 



^abs 



X e thr, 



(64) 



where we used equations ([H]) and (^) of Appendix. This ex- 
pression can be further averaged over the spectrum of scattered 
photons. The averaged values can be written as e±/c w p± = 

P±/n+. 
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FIG. 3. — Comparison of the analytical model (dotted curves) with the exact 
numerical results for a.2 = 2 and e^r = 10 2 . Here §ioad ~ 30 (see eq. W9I) 
and (acc = 5^ioad (see eq. pa). The break in n+ is at £br ~ 5 X IIP. 
The boundary of the ion-free zone § c is shown by solid vertical line. In this 
example, § c = 4 X 10 3 is chosen close to its typical value (see § 7). 

In the non-relativistic zone of the front the created e ± have 
initial Lorentz factors 7; n j = 1 + e±/2m e c 2 . Upon collec- 
tivization they push the medium forward and heat it. Using 
equations (^2j) and ( j43| ) we have 



7inj 



(1 + a^) 5 / 3 (a 2 



Ot\) 



2h+m e c (ct2 — on — 1)ckn 



(65) 



For hard spectra with finite ebr this equation overestimates 71 



(see discussion after eq. |5JJ); e.g. for a 2 = 1-5 and = 10 
the actual 7i n j w 9 (Fig. 1) while equation (65) gives 7j n j « 19. 
In the relativistic zone, £ ^> £ acc , we have 



X^thr 1 ( 4>w 



2m e c 2 7 27 7 V Ai 



l/a 2 



7 



< 1, (66) 



i.e. the Lorentz factor of created pairs is smaller than that of 
the medium. It implies that the e loading tries to decelerate 
the medium. Same effect can be viewed from the medium rest 
frame (see eq. p8| below). 
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To find the Lorentz factor of created e in the medium rest Suppose t r reaches unity at some w CT . Radiation scatters here 



frame, we use equations ( |l"T| ) and ( JT^ ) of Appendix and substi 
tute fi = ft, e sc = (1 + [3i)-\ and e thr = 2(1 + ft)/(l - ft) 
(cf. §§ 5.1 and 5.3). This yields 



7inj 



2 

2 



o(l- 

7 

Cthr 



/?ft)+Xethr(l 

, Xfthr _ _7_ 

4 7 ethr 



0)] 



» 1, 



(67) 



(we used 7 ^> 71). Here we keep the e sc term and neglect the 
e a bs term: in the rest frame, the scattered photon [blueshifted as 
7(27^) _1 > 1] is more energetic than the primary collimated 
photon [redshifted as (27)^]. In a similar way we evaluate 
momentum per injected particle in the medium rest frame, 



e thr> 

-7/ethr, 



(68) 



The e ± loading assists the medium acceleration only as long 
as pi n j > 0. At £ <; £ acc , p m j changes sign. At this point, 
ethr/7 ~ 1 and 7inj ~ (ethr/7) + (7/ethr) reaches its mini- 
mum ~ 2. The minimum of 7i„j is observed in Figure 1. 

Note that collectivization of the relatively slow e loaded 
at £ 3> £acc implies their fast acceleration (owing to e.g. 
stream instability, with the e ± stream being directed back- 
ward in the medium rest frame). In the case of partial collec- 
tivization (cf. § 2.2) the isotropized e acquire Lorentz factor 
7 e ~ 2jm e c 2 / e± 3> 1 found from equation (^) and then they 
are cooled by radiation. The cooling is dominated by photons 
with e ~ 1 and occurs in Thomson regime. When viewed from 
the lab frame, the hot move forward with Lorentz factors 
between 7 min = e±/2m e c 2 and 7 max = 2j 2 m P c 2 /e+. They 
scatter photons with energy e sc w e/2 ~ 1/2 (eq. [37]) through 
angles 1 — /z ^ (e±/m e c 2 )~ 2 — (x e thr)~ 2 - If the scattered ra- 
diation produced by this cooling process could interact with the 
primary beam we would have secondary pair production. This 
however requires 2/[e sc (l — //)] < eb r which in turn requires 
ethr < v^br/ f° r the primary injection event. This condition 
is not satisfied at ebr < 10 3 and the secondary pair production 
does not occur. 

6. BACKREACTION ON GRB 
6.1. Scattering 

The observed GRB can be affected by scattering in the cir- 
cumburst medium if the Thomson optical depth of the medium 
is substantial. In all calculations we assumed that the medium 
remains optically thin even after e loading. We now address 
this assumption. 

Consider a radius R and let the ambient medium have an ini- 
tial optical depth tr — no(R)aiR at this radius. In a medium 
with constant electron density no (R) = const, 



rji-7x 10- iu i?i 5 n . 
In a steady wind with mass loss M and velocity w, 



TR = 3xlO-^R-}^L. 

u e w 8 



(69) 



(70) 



The pair loading in the radiation front increases tr. The optical 
depth seen by the GRB photons at given zu is 



tr- 



n*(w) 
n 



(71) 



off the medium with 7„ = 7(w7 cr ) and acquires a new colli- 
mation angle 8 ~ 7~ 1 . This decollimation is not crucial if 7 cr 
is sufficiently large, 7„ > 7 m j n w 10 2 . Using the solution for 
n* /no and 7 (Fig. 1 and 3) one finds that the condition t r < 1 
at 7 = 100 reads 



tr < t c 



10- 



(72) 



This constraint slightly changes with a different choice of 7 m i n . 

Even assuming that the burst is sufficiently short (so that it 
decouples promptly from the ejecta, cf. case 2 in § 7.2) the 
whole burst can overtake the medium it sweeps only at R = R c 
(eq. ^3|). The condition (|72|) is therefore required at R > R c . In 
the ISM case (no = const ~ 1 cm -3 ) this condition is satisfied 
for any reasonable parameters. In the wind case the condition 
tr{R c ) < T cr reads (using eq. p3b 



M 21 . 7p 3/7 1/7 



A 

c 



(73) 



6.2. 7 — 7 absorption and the high-energy break 



In previous sections we assumed a priori a high-energy break 
in the primary radiation spectrum at some ebr- In a self- 
consistent situation, ebr is determined by 7 — 7 absorption of 
the primary 7-rays by the scattered radiation field. One can 
evaluate eb r in a simple way. At given tu the primary pho- 
tons e ~ e t hr are absorbed with rate h + (w). The number of 
absorbed photons during time R/c cannot exceed their initial 
number. This condition reads h+R/c < (Fi/m e c 3 )e^ h " 2_1 
and gives an upper limit on e tri r i-e. the self-consistent ebr- Us- 
ing equations (|58|), ([56]), and e t hr = 2(1 + ft)/ (1 - ft), we get 
(omitting a numerical factor ~ 1) 



ethr < 



1 



£1 n\Ra T 



£acc n-o 1 
£1 n\ tr ' 



(74) 



Hence eb r > 1 if TR(n* cc /no) < 1, i.e. the main radiation 
e ~ 1 is not self-absorbed after scattering if 



tr < 



n 



10" 



(75) 



This condition is weaker than the transparency condition ( \J2f) . 
To find ebr at 17? < 10~ 2 we need to solve the inequality ([74]) 
which is implicit since e t hr is a function of £1. The solution 
gives the maximum £™ ax and the corresponding e"^ x = eb r . At 
tr < 10" 2 we have e br > 1; then e thr w 87? = 8(£i/£ acc ) 6 
and n\ = 74/i e no(£i/£ aC c) 2 (using eqs. |62| and ^3]). We thus 
fmd£r7£acc = (600u e r Jt )- 1 /9 and 



0.1( Ue T R )- 2 ^, TR « 10" 



(76) 



The break appears in the GRB if zubi < A where zub T is 
given by equation (^0b. We have from equation ( |60| ) (using 
^acc ~ 120 A and eqT]l|) 



tU br W 6 X 10 



5:5 



6 x 10 S R 2 15 L 



53 



40 Q =(M e r fl ) 2Q 2/3 



(77) 
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The condition ro^r < A reads 



gbr 



2 x 10- 2 i?f,K 



15^53 



40^(M e ^i?) 2Q2/3 



< 1. 



(78) 



If the ambient medium is ISM with the optical depth ( |6^ ) 
then the condition (|78|) is not satisfied outside the emission ra- 
dius of the GRB and hence the GRB -medium interaction does 
not produce any break in the GRB spectrum. 

If the ambient medium is a wind with optical depth ([70]) then 
the condition (|78h is satisfied at radii R < R~ n where 



10 



15 



5.6 a2 u e 



3/(6+2a 2 ) 



Mai 
w 8 



2a 2 /(6+2a 2 ) 



(79) 



10 15 S 5 3 3 /10 (M 21 /^ 8 ) 2 / 5 cm. 

and cause a 



For instance, with = 2, i? 77 
It can be well outside the emission radius R c 
break in the GRB spectrum. 

We now derive e^r expected in the massive progenitor sce- 
nario. We substitute tr from equation ( f70| ) into equations ( f76| ) 
and ( fn\ ) and get 



<"br 



22R[ 



2/3 / ^8 



Af 2 



2/3 



(80) 



7? 



(6+2a 2 )/3 



n7 br w 6 x 10 s x 5.6 Q2 ^ 



^53 



\M21 



2q 2 /3 



cm. (81) 



The lowest ebr is produced at small R in the leading portion of 
the radiation front zu^t <§C A. With increasing R, vu^ and ebr 
grow. At R = i? 77 the tub r reaches A and ends up outside the 
front, and hence at R > i? 77 the 7-radiation does not change. 

From equations (80) and ( pl| ) one gets ebr oc vo^ +3 \ One 



can show that during the spectrum shaping at R < i? 77 the 
dependence ebr oc w 1 ^ 012 ^ is established at w < td^. At 
njbr < n7 < A a break is temporarily formed at a higher en- 
ergy ebr oc n7 1 /^ Q2+7 / 6 ^ (one finds the "extinction" zone 137 of 
primary photons e by equating their number to the number of 
scattered photons that have ethr = e)- At R > i? 77 , the whole 
front has e br oc w 1/( - a2+3 \ 

A distant observer will see first the leading portion and then 
deeper layers of the 7-ray pulse (the observer time is < Q b s = 
m/c). Eliminating R from equations ( |80| ) and (^Tj) we find 



£br(iobs) = (30 ± 10) (t ohs L 53 ) 



l/(a 2 +3) / W 8 



Al 



21 



2/(a 2 +3) 



(82) 



Here 30 corresponds to a 2 — 1-5, and the upper and lower val- 
ues correspond to ao = 1 and 2.5, respectively. 

The equation ( |82| ) has been derived assuming an idealized 
GRB with F{w) = const, while the observed GRBs are highly 
variable. The GRB variability did not play a big role to the 
structure of the front (see footnote in § 4); in contrast, eb r (^7) 
is sensitive to variability. The scattered radiation that shapes 
€b T (w) comes from a small leading portion of the GRB fluence. 
Therefore ebr is sensitive to the ratio of the flux at given w to 
that in the leading portion. The higher F{w) the less depleted 
its high-energy part is, which results in a higher e^iza). Hence 
a strong positive correlation between the instantaneous F and 
ebr should be observed in variable GRBs. 



7. EVOLUTION OF THE RADIATION FRONT 

In this section we study the front evolution with radius. We 
assume the internal scenario of GRB production (§1) and start 
from radius R cm where the emitted 7-ray pulse begins to de- 
couple from the ejecta. It happens when the ejecta becomes 
transparent, 



Rv.m — KR* 



6 x UpKt^E^T^-* cm. 



(83) 



Here i?» is the radius of "barion" transparency and K > 1 de- 
scribes a possible increase of the transparency radius owing to 
pair creation inside the ejecta; E c - 3 = r o jA/ C jC 2 is the energy 
of the ejecta. Note that K ^> 1 would require: (1) a substan- 
tial fraction of the emitted energy is in very hard 7-rays (above 
the threshold for pair creation, e > r o j), and (2) the hard 7- 
rays are emitted at a high rate at radii R^$> R* (otherwise 
production stops, pairs immediately annihilate to optical depth 
~ 1 and the ejecta becomes transparent on time-scale R/c be- 
cause of side expansion). It is unclear whether the two con- 
ditions are satisfied. The observed strong variations in many 
GRBs on time-scales St < 0.1 s suggest that in many cases 
R cm < 2T 2 c] c5t = 6 x 10 13 r e j2(^/0.1) cm. 

The thickness of the radiation pulse is equal to that of the 
ejecta, A rj A e j. Radiation is initially collimated within angle 
9 = r^ 1 and moves inside the ejecta. At R > R cm the radi- 
ation gets more collimated, 9 = r^ 1 (i?/i? om ) _1 , and gradu- 
ally overtakes the ejecta with relative velocity w (1 — /3 e j)c. 
The thickness of the radiation pulse emerging ahead of the 
ejecta and interacting with the ambient medium is growing, 
Ai(R) = c(l - Ay)(-R/e) = R/2Tl y When A t (R) reaches 
A e j, the whole pulse has left the ejecta. The corresponding £- 
coordinate of the back boundary of the interacting pulse £a is 



f « 1.1 x io A R^L 53 r e -- 

%i « 6.5 x 10 3 R- 5 2 E 53 , 



R 
2T5 

R 
2T 2 " 



> Aej- 



(84) 



7.1. i? cm < R < R S a±- Saturated surfing 

The pulse-medium interaction starts at R ^ R em with a very 
high £a ~ 10 6 . The medium entering the pulse accelerates to 
the equilibrium Lorentz factor 

7 sat (i?) = r ej ^- (85) 

at £ c ~ 10 3 <C £a and surfs the pulse. The acceleration time is 
~ (£ c A/c)7g at < R/c. Note that primary radiation is mixed in 
the front on scale <5n7 m ; x ~ 7sit^ because the photons have a 
finite angular dispersion 9 ~ 7^ . The Lagrangian coordinate 
tv is well defined only on scales 8w > zu m i X (on such scales 
the radiation can be assumed perfectly collimated with radial 
velocity c). Therefore, the ^-location of the medium in the front 
is defined with uncertainty £ m j x = n7 ln j X /A ~ £,a(R/ Rcm)~ 2 
which exceeds £ c at small R where £ c /£a < (Rem/R) 2 - 

The equilibrium Lorentz factor 7 sa t(-R) grows with radius. 
Correspondingly £ c [the value of £ where 7^) reaches 7 sat , 
see Fig. 1] grows and reaches ~ 10 4 at R ~ 10 14 cm. At 

R = RsaU 



R sat w 1.3 x 10 u T ei - 3/i L^ReJ/* cm, 



(86) 



7 sa t exceeds 7 max given by equation (^9)). Then the medium 
cannot accelerate to 7 sa t on time R/c and the saturated stage 
ends. 
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7.2. R sat < R < -Rgap- Unsaturated surfing: caustic 

Now £ c and -f c are determined by the condition (A£ c /c)7^ ~ 
i?/c (the time of acceleration to -f c is about i?/c). Using equa- 
tion (p2|) we get at 7 C > 27, 



53 z- 1 / 4 ;?- 1 ^ 

1/4^53 ^1 
4"acc 



l 15 



7c 



2.0 x 10 3 

TVs 

Sacc 



r 3/8 p -3/8 
-^53 -"15 



(87) 

(88) 



When R grows from 10 14 cm to 10 6 cm, £ c decreases slowly 
from £ c ps 30^ acc to £ c ps 10£ acc . Correspondingly, 7 C de- 
creases from ps 10 3 to ps 140. Hereafter we substitute in all 
estimates £ aC c = 120 keeping in mind the typical «2 = 1.5; for 
oi2 = 2 there is a slight change: £ acc = 150 (§ 5.2). 

The new material trapped at given R comes to £ c with smaller 
7 compared to that of the already accumulated material in the 
front. This results in "overshooting" and implies appearance 
of a caustic. The overshooting can be seen e.g. in the in- 
coordinate: the accumulated material has ti7° ld = £ C A and 
the newly trapped material comes to nv" cw = R/j c - Hence, 
va^ cw /va° ld oc £c/7c With decreasing £ c , the condition for 
overshooting reads d(£/7) /d£ < at £ = £ c , or 



din 7 
dln£ 



> 1, 



(89) 



which is satisfied (see eq. |62|). The caustic results in a shock. If 
the shock is radiative (which may be the case since the material 
is pair-dominated and the Compton cooling is very efficient) 
then the shocked matter piles up in a thin shell. 

When the caustic appears, the accumulated ion material be- 
gins to decelerate and the stream loaded behind £ c hits the 
ion medium. One can show that the momentum of the 
stream exceeds the momentum of the accumulated ions and this 
"reverse" shock should be strong. For simplicity, we will ne- 
glect the impact of the e ± stream on the surfing medium (the 
inclusion of this effect will slightly increase the radius i? gap 
derived below). 

The medium surfs the pulse with 7 ps 7 C until £a reaches £ c . 
This happens at some radius R c . We now evaluate R c in two 
possible cases. 

1. R c < 2r^jA c j. — The ejecta catches up with the surfing 
medium before the whole 7 -ray pulse leaves the ejecta. Then 
£ c = £a gives (using eqs. p4jandl 



Rc 



1.2 x 10 



18 



A") 3 



8/3 
ej2 



cm 



10 1 



-53 



8/3 

cj2 



cm. 



r 



(90) 



At R = R c we also have ~f c ps T C j i.e. the ejecta touches the 
medium with a small relative Lorentz factor and starts to de- 
celerate. The gap between the surfing medium and the ejecta 
disappears at this moment; we thus have i? gap = R c 
The assumed condition R r 



< 2r^A cj requires t b > t%, 



45 E, 



1/2 
53 



A/2 p 

Sacc -L 



7/3 
oj2 



4 ^53 



1/2 



r 



7/3 



(91) 



<\i 2 



where we used tf, = A C j/c = E/L. This condition implies 



-Rgap Rc 



< 



3 x 10 16 E, 



1/2 
53 



A/2 

Sacc 



1/3 



3 x 10 15 — ^- cm. 



CJ2 



1/3 



C J2 



(92) 



2. R c > 2r 2 .A ej . 



- The whole 7-ray pulse leaves the ejecta 
before they reach R c . Then £ c = £a at 



Rc 



1.6 x 10 1 

Sacc 



17.3/7,1/7 



2 x 10 15 E^ 7 ^ /7 cm. 



The value of 7 c (i? c ) now differs from r c j, 



7c 



7.1 X 10 2 p3/ U ,-3/7 
_£/ 53 Z >- 



>3/14 
Sacc 



(93) 



(94) 



The condition i? c > 2r 2 jA C j (which is equivalent to t\, < t%, 
cf. eq. [h]) implies that -f c > T c j and hence the gap still exists 
at R = R c . The gap disappears only when 7(£a) ~ r c j (then 
the ejecta catches up with the surfing medium). This condition 
yields 



Ra 



3.0 x 10 16 E, 



1/2 
53 



A/2 

Sacc 



r. 



1/3 

cj2 



cm ps 2.7 x 10 1 



E 



1/2 
53 



1/3 



cm. 



CJ2 



(95) 



> tt) and case 2 



Note that R gap is smaller in 



We call case 1 as "long-burst" regime (tb 
as "short-burst" regime (tj < ti) 

case 1 (compare eqs. ^5| and ^2|). The gap may be not opened 
if the GRB emission radius is much increased by e ± inside the 
ejecta (cf. eq. |83|) and exceeds i? gap ; -Rem > -Rgap would re- 
quire K > 170r ej2 " 2/3 (E/E ei ). 

Hereafter in this paper we focus on the short-burst regime. 
Then the energy of the 7-ray pulse interacting with the am- 
bient medium is constant at R > R c and the simple scaling 
£a oc R~ 2 holds. A simplified blast wave model can be con- 
structed for short bursts (as done in § 8). The extension to long 
GRBs is straightforward, though it implies additional technical 
details which we defer to a future paper. 

7.3. i? gap < R < i? acc - Preaccelerated pair-rich medium 

In this range of radii, £a decreases from £c(-Rgap) ~ 10 3 
to ^acc ps 10 2 . Correspondingly, 7(£a) decreases from r e j to 
ps 1. When £a < £acc the front cannot accelerate the medium 
to relativistic velocities. The condition £a = £acc therefore de- 
fines the maximum radius where the relativistic preacceleration 
occurs, 



i? acc - « 7 x lO 15 ^ 2 cm. 

4~acc 



(96) 



7.4. i?acc < R < i?ioad- Non-relativistic pair-rich medium 

At R > R acc the radiation front still loads the medium with 
a large number of pairs. At R = i? acc (£a = £acc) we have 
n* /no ?s 74/i e behind the front, and with increasing R the pair 
loading decreases exponentially (see eq. [|9]). The pair loading 
ends at R = R\ oa a (here £a reaches £ioad and n* /no ~ 1), 



-Rload — 



Rx 

A/2 
Moad 



1.6 x lO 16 ^ 2 cm. 



(97) 



In § 5 we showed that £ acc is related to ^load by a simple for- 
mula £ acc = (5 + ln/i e )£ioad which weakly depends on \i e 
(1 < p, e < 2). Hence, we have a relation 



i?load = (5 + m/Xe) 1/2 #acc = (2.3 ± 0.1)i? a 



(98) 
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7.5. R > -Rioad- Front weakly affects the medium 

Here £a < £ioad and both loading and preacceleration are 
shut down. The blast wave sweeps the normal pair-free medium 
which has 7 w 1. 

8. BLAST WAVE 
8.1. Dynamics 

When the gap is closed, the ejecta drives a blast wave through 
the medium preaccelerated by the leading radiation front. We 
will model the blast wave in a simplified way, as a thin shell 
sweeping the medium. This is a good approximation to the ex- 
act hydrodynamic solution with forward and reverse shocks if 
the ejected shell is sufficiently thin, so that the reverse shock 
crosses A e j on time less than R/c (e.g. Piran 1999) 4 . 

The shell has initial mass M R \ and Lorentz factor T e j and 



= -Rgap (§7). At a radius 
> M c \ and Lorentz factor 



starts to sweep the medium at R 
R > i? gap the shell has mass M 
r < r c j . When it sweeps a mass element dm that moves with 
Lorentz factor 7, T decreases by dT and energy di?diss is dissi- 
pated. The laws of energy and momentum conservation read 

YM + 7dm = (r + dT) (M + dm + dm hcat ) , (99) 
T0M + j0dm = (T + dT)(0 + d/3) 

x (M + dm + dm hcat ). (100) 

Here = (1 — 1/r 2 ) 1 / 2 is the shell velocity and dmhcat = 
(di?diss/c 2 r) is the rest mass associated with dissipated heat. 
The inertial mass M includes the initial mass of the ejecta 
Moj, the swept mass m(R), and the stored heat. We will as- 
sume that a fraction 77 of dmheat is radiated away and the 
rest remains to increase the kinetic energy of the shell. Then 
dM = dm + (1 — ?/)dTOhoat and we get dynamic equations, 



M- 



dM 
dm 

The radiated energy is 
c 

'/- 



dT 
dm 

(1 



= r 2 /3 7 (/3-/3), (101) 
77)r 7 (l -00). (102) 



dm 



dm 



= T]C Z T 



T~/(l-00)-l . (103) 



The swept mass is related to radius by dm/dR = AttR 2 po 
where po(R) is the medium density ahead of the radiation front. 
Note that the increase of the medium mass by e ± loading is 
fairly low at R > i? gap (§ 7). The front affects the blast wave 
dynamics mainly by increasing 7 of the ambient medium (and 
possibly by increasing the efficiency 77 as a result of e ± load- 
ing, see TM; here we assume r\ — const for simplicity). The 
dynamic equations acquire the standard form if 7 = 1 (deceler- 
ation by static medium, see Piran 1999). 

The time interval between radiative preacceleration and sub- 
sequent sweeping by the ejecta shell is much smaller than R/c, 
and hence preacceleration can be treated locally at a given R. 
Indeed, the distance between the leading boundary of the radi- 



ation front and the blast is A 
sweeping time is t s 



f = (1 - 0)R w R/2T 2 and the 
(A / /c)2 7 2 = (i?/c)(7/r) 2 < J R/ C at 



pulse interacts with the ambient medium ahead of the ejecta (we 
assume the short-burst regime, so that the whole pulse has left 
the ejecta at R < R gap , cf. § 7.2) and 



Rx 
R 



U = [ — ) = -5- 



R 



We will use the analytical formula 
1, 



Rs 



for 7(£a); then 



(104) 



7(x) 



-1/3 „ 



X > 1, 

71 <a;<1 ' 



(105) 



< X < 



1 

71' 



where x gap = y^ST^ 1 ''' w 0.3 is found from 7 = r oj . 

The characteristic mass of the problem is the ambient mass 
within the acceleration radius, 



4irR 2 p dR. 



(106) 



The mass swept when the blast wave reaches a radius x is 



m{x) 



mjvrrX 



(107) 



where k = 3 for a constant-density medium and k = 1 for a 
wind with constant M and w. 



At 7 <C T equation (103) yields 



dm 



s 2 P 



(1 + 0) 



7 



(108) 



TciMcjC 



Before the ejecta decele rates , i?diss <C E c 
have r w r o j. Equation (108) then yields [we use eq. |105| for 7 
and replace 1 + by the step function 1 + H (1 — x)] 



-EdissOr) 



T^2 

r ej m a 



x - 

2kx 6 



6-k 
6+fe 



6+fc 
2k(x 3+k - 



+ "0, X > 1, 

•0, 



< 1, 



' gap 



< X < 



1 

73' 

(109) 



where t/j « 0.004 for k = 1 and i\) w 0.002 for k = 3. Equa- 
tion (|109|) assumes a deceleration radius Xdcc > L Setting 



E e ; we find the actual deceleration radius, 



•^dc 




D 



l/k 



2^r 2 jm acc c 2 
(Q + k)E ci ' 



D < 1, 



1 < D < 27(v / 3)* 



(HO) 



(HI) 



Note that Xdcc depends very weakly on D at D > 1. Equa- 
tion (|l05l) gives 7(x doc ), 



7dc 



1, 

D d/(d+k)_ 



D < 1, 

1< L> < 27( v / 3) fe . 



(112) 



i? > -Rgap- The Lorentz factor 7 at given i? is that found behind 
the radiation front, 7 = 7(£a) (§§ 4 and 5). The whole 7-ray 

4 The standard blast wave model assumes the formation of collisionless shocks. Smolsky & Usov (2000) developed and alternative model for the ejecta-medium 
interaction. The sweeping-shell approximation is useful in that case as well as it deals with energy-momentum conservation only and gives a correct -Rdcc- 



The simple formula ( 1 1 C ) neglects ip and applies if D < 
27 (v / 3) fc (corresponding to Xdcc > 1/VS). The extension to 
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D > 27(v / 3) fc is straightforward. A much higher D, however, 
would imply that the swept medium is optically thick, making 
the model inconsistent (we assumed transparency when build- 
ing the front model in §§ 2 and 3). Let us evaluate the optical 
depth of the swept medium at Rgag- The e ± loading factor is 
n*/n a pa 2 x 10 2 /j e /cc 2 ap (cf. eq. [53|), and 



2 X 10 2 TOgapCT 2 X 10 2 TO acc <7T 



x gap 47ri? 2 ap m p 



4-fc 
(-gap 



47ri? 2 cc m p ' 



(113) 



The condition t* _ < 1 defines a maximum Substituting 

gap ■ © 



i' gap 



= \/3r • and using equations (pq) and (1 1 1 ) we find 



D. « 0.1 



cj ^ -p(2+fe)/3 



(114) 



With a typical r cj = 300, pa 30 for k = 1 and D* pa 10 3 
for = 3. If £) > £>* the trailing photons of the radiation front 
will be trapped and advected by the medium until it becomes 
optically thin. Our front model does not apply to this case. 

It is in stru ctive to compute the swept mass at i?dcc (using 
eqs. 11C and 107). In the regime D < 1 {xdec > 1, 7doc = 1) 
the standard estimate holds (e.g. Rees & Meszaros 1992), 



m dc 



D < 1. 



(115) 



In the regime D > 1 (xdec < 1, 7dcc > 1) the deceleration 
occurs in a relativistically moving medium; then 



TOde 



(6 + k)E cj 6 /(6+k) = ( 1 

2fcr 2 c 2 ^ 



6 + k \ 7doc^e, 



2k 



If the ambient medium is ISM with a constant density no 



1— 10 2 cm 3 we have k = 3 and m acc 
The parameter D is then given by 



= 2.4xl0^ e «o^53 2 g 



D w 1.4 x 10 



~ 4 fj, e n E^ 2 E, 



l-p 2 
6J53 e J2- 



(117) 



The t ypica l D < 1 and hence Xdcc > 1 and the standard esti 



mate (115) applies. Equation (109) yields the energy fraction 
that is dissipated at x < 1: / acc ~ (2/3)x^ ps -D < 1. 
The fraction dissipated in the static pair-loaded zone, 1 < x < 
zioad « v 7 ^ (eq. ||), is / load « (xi oad /a;dcc) 3 « 20D. 

If GRBs have massive progenitors (e.g. Woosley 1993) 
their ambient medium is a wind from the progenitor. From 
a Wolf-Rayet progenitor one expects a wind with mass loss 
M — 10~ 5 M Q yr" 1 and velocity w — 10 3 km s -1 (Cheva- 
lier & Li 1999). In the case of a red giant, the wind veloc- 
ity is smaller, w ~ 10 km s _1 , and then the ambient den- 
sity is higher. The wind medium is described by k = 1 and 

m aC c = {M/w)R acc = 7 x 10 28 M 2 iu>g x £^ 2 
D is then comparable to or much larger than unity, 



2iw s 'iigg g. The typical 



D = 2Mc 2 r c / 



7wE c 



l.%M 21 w^Eli 2 E c ^Y c ^. 



(118) 



Besides, a s mall l ocal maxi- 

" (l + /?)/7 



108 



dissipation rate peaks at x pa Xd. 
mum appears at x = 1 [it is understood from eq 
has a maximum at j3 = 0.5 i.e. at x = 1]. If D > 1 then Xd cc 
is close to unity and one can see the strong peak of energy dis- 
sipation at x ~ 1. For Z? = 100, 80% of the blast wave energy 
is dissipated at 0.5 < x < 1 and 99 % at 0.3 < x < 2. Note 
that in this model D > pa 20 (eq. 1 14) and the initial stage 
0.3 < x < 0.5 is optically thick. Yet, since most of the energy 
is dissipated when the medium becomes transparent, the model 
with D = 100 is marginally applicable. 
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FIG. 4. — Blast wave dynamics in a wind. Here E = E e j = 10 53 erg 
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The dissipated energy Z?diss (all emitted if 77 = 1) exceeds 
E e - 3 because the additional energy _E acc = ° 2 J(l~ l)dm is de- 
posited by the radiation fron t whe n it preaccelerates the ambi- 
ent medium. Using equation (105) one finds for a wind medium 

(k = 1) 



E n 



' 3^3 

2a; gap 



5.9 




5.9 



(119) 



To study the wind case in more detail we solve numerically 
equations (101,102,105,107). Figure 4 shows the results for 



r cj = 200, E = E cj = lO^erg, and the efficiency 77 = 1. The 



The main contributio n to E acc comes from small radii x ~ 
£gap- Using equation (|1 1 1|) one can find E acc /E e j as a function 



of D. For example E a 



2.1 x 10- 3 £>forr cj = 200. 
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The parameters r\ = 1 and E c - } = E taken in the numerical 
examples imply that the blast emits energy equal to that of the 
prompt GRB. This emission should contribute to the preaccel- 
eration (it is soft and preaccelerates efficiently, with no Klein- 
Nishina reduction). Here it was not accounted for; in the model 
with D = 100 it would increase i? acc by a factor of 2. 

8.2. Emission 

We now evaluate the main characteristics of the blast wave 
emission, in particular, the bolometric light curve seen by a dis- 
tant observer and the synchrotron peak frequency. A detailed 
analysis is deferred to a next paper. 

8.2. 1 . Bolometric light curve 

A distant observer will see a mixture of radiation emitted by 
the decelerating shell at different radii. Denote the arrival time 
of radiation by < Q b s and choose t ^ s = for a light signal that 
would come from the center/beginning of the explosion. First 
consider the observed light curve from instantaneous emission 
of energy Eq by the shell at radius R. The shell reaches this ra- 
dius at a time t(R) after the beginning of the explosion. Choose 
a polar axis 9 = pointing toward the observer. The observer 
will first receive photons emitted at 8 = (// = cos 9 = 1). 
These first photons come at i b s = t(R) — R/c and photons 
emitted from a circle fi = const < 1 arrive with a delay of 
(R/c)(l — p). We thus have a relation 



R 

t obB (-R,M)=*(ir)--M, t(R) 
c 



dR 

o (3c 



(120) 



Radiation received in time interval <5i bs comes from the ring 
\5p\ — (5t bs(c/i?). The total energy emitted by this ring (in all 
directions) equals SE = Eo\5p\/2 where \5p,\/2 is the frac- 
tion of the shell surface occupied by the ring. We will as- 
sume that each element of the ring emits isotropically in its rest 
frame. Radiation emitted toward the observer within a rest- 
frame solid angle dSl occupies dSl = T 2 (l — (3p) 2 dVl in the 
lab frame. Hence the observed flux is affected by the beaming 
factor r~ 2 (l — (3p)~ 2 and the apparent isotropic energy seen 
by the observer from a ring Sp is SE app = T~ 2 (l — j3p)~ 2 8E. 
(Integration of SE app over the shell gives Eq as it should be.) 
The apparent isotropic luminosity is 



_ SE app 



St. 



obs 



(121) 



From the dynamic solution we know dE^iss/di? and T(R). 
It allows us to compute the observed light curve from the 
whole history of the shel l de celeration [we substitute Eq = 
r)(dEdiss I dR)dR into eq. 121 and integrate over R], 



L bs{t bs) — 



{R/c)r,(dE diss /dR)dR 
2T 2 [R/c-^(t-t ohs )] 2 



(122) 



where R m „^ (t n h s ) is defined by condition t(R) — R/c = t ^ s 



(see eq. |120[ ). 

The results are shown in Figure 5 for the blast waves in wind 
environment (same cases as in Fig. 4). There is no emission 

until t risc = i? gap /2r 2 j a E^ 2 (T C -J2QQ)- 7 ^ s which cor- 
responds to the moment when the ejecta catches up with the 



surfing medium and begins to decelerate. At t ^ a 
light curve rises steeply and then reaches a peak at 



''peak 



2^ 



1- 

■^dec ) 



D 
D 



< 1 
> 1 



1/2 



CJ2 



irise the 



(123) 



Since Xd cc remains close to unity even at D ^> 1, we get a 
universal t pca k in a very wide range of D. In the examples 
shown in Figure 5 (E 53 = 1, r e j 2 = 2) we get i pca k ~ 3 s. At 
D < 10~ 2 there appears a plateau in the light curve between 
tp Ca k and ~ 0.3i?dcc/2r 2 jC. In the model with D = 100 the 
total emitted energy exceeds E C j by 20% owing to E acc (§ 8.1). 
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FIG. 5. — Bolometric light curves from the blast waves shown in Figure 4. 

Figure 5 assumes a radiative blast wave (rj = 1). Adiabatic 
blast waves (77 <C 1) produce similar light curves with the same 
ipeak- In reality the radiative efficiency is likely to change dur- 
ing the afterglow, being highest at R ^ R acc (when a lot of 
are loaded) and decreasing at larger R. Then the peak will be 
sharper. 

The light curve should further be corrected for the cosmolog- 
ical effects. In particular, for a burst with a redshift z one gets 

Vak » 12 (! + ^53 /2r W 2 S - 

8.2.2. Synchrotron peak frequency 

The medium encountered by the blast wave is e ± -loaded, 
preaccelerated, and compressed by the leading radiation front. 
Correspondingly, the standard analysis of the blast wave emis- 
sion (Blandford & McKee 1977; Piran 1999) applies to our case 
with three modifications: (1) er- loading increases the number 
of shocked electrons by the factor n* /no. (2) The proper mass 
density of the medium is p = /9o7 _1 (l — (3) rather than pq\ 
here the factor (1 — 0) comes from compression in the front 
(cf. § 3.2) and 7 _1 appears due to Lorentz stretching when we 
go to the rest frame. (3) The Lorentz factor of the pre-shock 
medium in the shock frame is r/7 rather than T. 

The proper energy density of the post-shock material is 



I I ~- ] pc 2 
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Assuming that magnetic energy B 2 / 
we have 



iir is a fraction £g of a, 



B 



£ / 32ire B po 



(125) 



Assuming that e share the energy of shocked ions, the mean 
randomized Lorentz factor of e ± in the rest-frame of shocked 
matter is 



7e 



m* r 

ra e 7' 



(126) 



where = u e m p no / n* and 1 < (i e < 2 (eqs. [13] and [l4])- 

We then get the peak synchrotron frequency in the rest-frame 
of the shell, 



lO b B~/ 2 Hz 



10 30 po /2 



7 



5/2 



^ e i/2 - 3 



ri Hz. 



(127) 



The corresponding observed frequency is v s = i> s r(l + z)^ 1 
where z is the cosmological redshift of the burst. 

For example consider a blast wave with D <C 1 in ISM. At 
radii i? acc < R < Rdec we have T r o j and 7 « 1. The 
density n* is given b y equ ation ( p^ ) with ur/a = {R\ oa d/ R) 2 - 
Then from equation ( 127 ) we get 



5 x 10 V /a ei /a nS /a ^» 



exp(2i?Ld/^ 2 ) + cx P (-2 J Rf oad /i? 2 ) + 2 

Race < R < Rd 



Hz, (128) 



i? acc has the 



For instance, the observed emission from R 

peak frequency v s w 2 x 10 16 (1 + z)~ 1 e 1 J 2 /i e 
The emission from i? < i? acc is even softer, however, the lu- 



acc 

1 1/2 1/2 l/2 r 4 



r ej2 Hz. 



D at i? acc , see 



minosity is small from that region (£^1 
§8.1). 

As a second example consider a blast wave in a wind with 
D > 1 and evaluate the peak frequency at the dece leratio n ra - 
dius i?dcc < Ra.cc - We substitute equations ( |1 10| ) and ( 112| ) 
into equation (127) and use n* /no ~ 74u e x~ 4 (cf. eq. p3); po 

w) and 



is expresse d in terms of D using po = (M/AttR 



equations (118) and (661). Then we get 



acc dec 



2 x 10 
D 2(i 
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For instance 

v s = u a T{\ 



with D = 10 and = 0.1 the observed 

-1 



zj "is close to optical/UV and then a large 
fraction of the blast wave energy is emitted in the optical band. 



8.3. Comparison with the standard afterglow model 

The impact of the radiation front on the blast wave is illus- 
trated in Figure 6. We take the usually assumed parameters of 
a GRB with a Wolf-Rayet progenitor and compare the model 
developed in this section with the standard model that neglects 
the radiation front and has 7=1 and no ahead of the blast 
wave. For a simple illustration we assume a short GRB with 
an impulsive ejection (so that the sweeping shell approxima- 
tion can be used, cf. § 8.1). The synchrotron peak frequency 
v s is evaluated with es = 0.1 which is close to its maximum 
possible value. 




FIG. 6. — Afterglow from a GRB ejecta decelerating in a progenitor (Wolf- 
Rayet) wind with M = 2 X 10 _5 M Q yr~\ w = 10 3 km s" 1 , and /i e = 2. 
The GRB is modeled as an impulsive emission of a gamma-ray front (with 
isotropic energy E = 10 53 erg) and a thin ejecta shell with kinetic energy 



10- 



imply D ; 



' erg ar 
10 (eq. 



4-L 



orentz factor Y, 



200. The assumed parameters 



118). Dashed curves show the prediction of the standard 
model that neglects the impact of the radiation front and solid curves show the 
actual behavior. Two extreme cases are displayed in the figure: r\ = (adi- 
abatic blast wave) and 77 = 1 (radiative blast wave). Four zones are marked: 
I — R < Rgap (the gap is opened), II — J? gap < R < J? acc (the gap 
is closed and the ejecta sweeps the relativistically preaccelerated er^ -loaded 
ambient medium), III — -Race < R < R\ oa ,d (e^-loaded ambient medium 
with 7 1), and IV — R, > Ri oa d (pair-free ambient medium with 7 ss 1). 



Radius is measured in units of i? a 



0.7 x lO 16 ^ 2 cm. Top panel: the 
dissipation rate. Bottom panel: the synchrotron peak frequency u s (assuming 
eg = 0.1) in units of m e c 2 /h. No cosmological redshift correction is applied 
here; the observed u s should be further reduced by(l + ^) _1 . 

The front impact is dramatic in both radiatively efficient 
(7/ = 1) and adiabatic (n -^i 1) cases: (1) no afterglow is emit- 
ted at R < Rgap, (2) the afterglow peaks sharply at R = R acc , 
and (3) the main emission occurs in the pair-loaded zones II and 
III and it is much softer than predicted by the standard model. 
Note that v s is highest at R w R acc for n = 1 and at R w i?i oa d 
for n — 0. In the former case v s lies in soft bands at all R. 
For smaller eb and larger D the afterglow will be even softer. 
One thus can expect "X-ray-weak" afterglows from GRBs with 
massive progenitors. 

Note that the basic dynamic equations of § 8.1 are not accu- 
rate in the adiabatic case. If not radiated, the dissipated energy 
is subject to adiabatic cooling and tries to accelerate the ejecta 
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(or rather to reduce the deceleration rate). The present formu- 
lation neglects this effect and is in error by a factor of a few. 
Another special feature of the adiabatic regime is that the time- 
integrated dissipation rate far exceeds E e] . The shell kinetic 
energy remains approximately constant and it is redissipated 
many times. 

9. CONCLUSIONS 

The unusual character of radiation fronts in GRBs is owing 
to two basic facts: (1) the front is opaque for ambient electrons, 
i.e. the electron scatters many times when it is overtaken by the 
front, and (2) the front is opaque for the decollimated scattered 
photons (7 — 7 opacity). The two properties cause e ± load- 
ing and violent acceleration of the medium. The important role 
of the front is that it "prepares" the medium encountered by the 
blast wave. We summarize the medium dynamics in the front in 
§9.1 and its impact on the ensuing blast wave in § 9.2. The front 
should cause spectacular observational effects during the early 
afterglow. The expected phenomena are discussed in § 9.3. 

9.1. The radiation front 

1 . Photons scattered at one portion of the front get absorbed 
at a different portion behind the location of scattering. 
The local approximation used previously (assuming that 
the scattered photons instantaneously become e ) is not 
adequate: the front structure is essentially governed by 
the radiative transfer. Yet a simple analytical description 
can be given to this non-local structure (§ 5). 

2. At radii R > R c (§ 7.2) the whole front has a quasi- 
steady structure established on time-scales <C R/c (an 
ambient particle spends time St <C R/c in the front). The 
propagating front is described by a self-similar solution: 
the density amplification nj no and the medium Lorentz 
factor 7 are unique functions of £ = wj A. Here < 
-uj < A measures the distance inside the front (w = at 
the leading boundary) and A oc R 2 is the electron free- 
path in the radiation field (eq. [j]). A simple analytical 
approximation to the structure functions is given in § 5. 
In the leading portion of the front the medium density ex- 
ponentiates due to pair loading on length £i oa d ~ 30, at 
£acc ~ 5£ioad the medium is accelerated relativistically, 
and at £± w 30£i oa d the loaded pairs outnumber the am- 
bient protons by the factor m p /m e and dominate the in- 
ertia of the medium. The medium parameters behind the 
front are «(£a) and 7(£a) where £a = A/A oc R~ 2 is 
its trailing boundary. The front acts as a relativistic ac- 
celerator at radii R < i? acc = 7 x 10 15 i?g3 2 cm where 
£,A > £acc- At R > i? acc the front still loads the medium 
with e ± . At R > i?i oad = 1.5 x 10 16 £^ 2 cm the e ± 
loading is shut down (£a < £ioad)- 

3. At radii R < R c the medium is accelerated so strongly 
that it gets "stuck" in the radiation front (the time-scale 
for the medium dynamics across A exceeds R/c). Then 
two zones exist in the front: (1) £ < £ c ~ 10 3 — 10 4 — 
here the steady self-similar structure is established and 
(2) £ > £ c — the ion-free zone. Being strongly accel- 
erated, the ambient medium cannot penetrate the zone 
£ c < £ < £a- Instead, it accumulates at £ ~ £ c and surfs 
the radiation pulse. With increasing R the front traps new 
material which is accelerated to a smaller velocity. This 
causes the overshooting effect and a caustic appears in 



the surfing medium (§ 7.2). The surfing stage ends at the 
radius R c which is equal to or smaller than the radius of 
the blast wave formation. 

9.2. The blast wave 
1. The blast wave forms at R = i? gap . For short bursts 



3 x io 15 i4 /2 r ej - 



1/3, 



(t b < AEli 2 Y Ci2 V \) i? gap 

and for longer bursts i? gap ~ 10 16 t b " 1 i?53r e j^ 8 ^ 3 cm. 
We limit our consideration here to short bursts. Then the 
7-ray front is detached from the blast wave (by a small 
distance R/2T 2 ) and it leaves behind the changed ambi- 
ent medium as described in § 9.1 (item 2). In contrast, 
radiation of long bursts continuously "leaks out" of the 
ejecta during the blast wave stage which implies addi- 
tional technicalities deferred to a future work. 

2. If the explosion happens in a constant-density medium 



Withpo < 2 X 10- 22 r ej2 2 £e J53 -53 



eJ 12 g cm" 3 then the 



blast wave decelerates at i?dec > -Rio a d where the impact 
of the radiation front is small and the standard afterglow 
model works well. Such conditions probably take place 
for explosions in ISM. 

If the explosion happens in a wind from a massive pro- 
genitor the radiation front affects strongly the b last wave 
deceleration. We defined a parameter D (eq. 118) that 
controls the dynamics; the blast waves in winds have 
D > 1 which corresponds to i?d cc < fia.cc i.e. the 
deceleration occurs in a relativistically moving medium. 
The standard estimate of Rdc. r i s then invalid and instead 



one should use equation (110). The ejecta decelerates 
near the unique radius i? acc w 7 x 10 15 i?53 2 cm; it im- 
plies that i?dcc weakly depends on the medium parame- 
ters as long as D > 1. The ejecta does not decelerate 
at R <C Race because the medium ahead has a high 7; 
at R i? acc the medium preacceleration ceases rather 
abruptly, 7 = (i?/i? acc )~ 6 , and the delayed deceleration 
occurs violently. 

The strong effect of the radiation front on the blast wave for- 
mation and dynamics is easily understood. The front passes en- 
ergy SE w StE/2 to an ambient mass Sm ahead of the ejecta, 
where St w 0.2<7t^to(47t i? 2 m*) _1 is the optical depth of Sm 
(here 0.2 is a Klein-Nishina correction and to* is mass per elec- 
tron, < m p due to pair loading). As a result the medium is 
accelerated to a Lorentz factor 7, 



7 -l = — = 
c 2 Sm 



St E 
Sm c 2 



O.Ictt E 

47ri? 2 m* c 2 



(130) 



The front structure solution gives to, and 7 behind the front; 
e.g. to* 1=3 TOp/74 at R = i? acc (/3 = 0.5). At small radii, 
R < R g n P ~ i?acc/3, one gets 7 > r o j, i.e. the accelerated 
medium runs away from the ejecta and the gap is opened. 

The ejecta decelerates when it sweeps a sufficient mass to. 
Namely, the deceleration condition reads (r e j/7)m w M e j i.e. 
the swept inertial mass measured in the ejecta frame equals 
M e j . The radius of deceleration is increased if 7 3> 1 despite 
the fact that only a small energy e -C E C j = T C jM C jC 2 was 
used to preaccelerate the ambient medium. Indeed, we have 
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[omitting the numerical factor (6 + k)/2k, cf. eq. 116], 



m dc 



7doc M _ 

1 pi 



(7dec - l)m de 



E r 



7doc(7dcc - 1) 



r 2 



(131) 



The growth of m^cc oc 7dcc occurs while e < E C j for 7d cc < 



-ej- 



9.3. Expected observational phenomena 
1. The generic prediction is that the early emission of a 



GRB blast wave (at t ohs < 30E, 



1/2 
53 



r - 2 



s) should be 



very soft 5 . Compared to the standard model that neglects 
the effects of the radiation front, the peak frequency of 
synchrotron emission is reduced by the pair loading fac- 
tor (m*//K e m p ) 2 = (n*/no)~ 2 and the preacceleration 
factor ^~ 5 / 2 (see eq. 127). The early afterglow should 



start as a relatively weak optical emission at R < i? acc 
and then the peak frequency moves to the X-ray band; 
only at R > i?i oa d the blast wave sweeps the normal 
e ± -free medium with 7 w 1 and emits in the standard 
regime. 

2. The fraction / of the afterglow energy that is emitted at 
the early soft stage is controlled by the ratio i?i oa d/-Rdcc- 
In the typical ISM environment / < 10%. In the typical 
wind environment (massive progenitor scenario) / ~ 1 
i.e. most of the blast wave energy is emitted at the early 
soft stage. In the latter case we emphasize the likely 
possibility of "X-ray-weak" afterglows whose emission 
peaks in soft bands throughout the whole evolution of 
the blast wave (Fig. 6). 

3. The light curves from blast waves in winds have spe- 
cial features which may be easily recognized in observa- 
tions. The violent deceleration that happens at R ~ -R aC c 
should cause a strong peak in the soft light curve. Then 
one should observe (1) a steep rise of the afterglow at 



tr 



-Rga P /2r 



1/2 2 

i? 5 3 r c j 2 and (2) a peak at 



2 - 1/,3 f 

ti.OL e j 2 l-r 



Both t r 



^pcak ~ {E acc / Rg a p^)t T [ st 

and < pca k depend weakly on the wind parameters in a 
wide range 1CT 3 < D < 10 2 (§ 8.2.1). Given the ob- 
served E and £ pca k one can find the Lorentz factor of the 
ejecta. If T e j does not vary strongly from burst to burst 
(as suggested by the clustering of GRB spectral peaks at 
e ~ 1, see Preece et al. 2000) there should exist a strong 
correlation between i poa k and the observed isotropic en- 
ergy E of the prompt GRB. 

4. In the massive progenitor scenario, the prompt high- 
energy 7-rays must be absorbed efficiently by radiation 
scattered in the wind (§ 6.2). As a result, the high-energy 



tail of the GRB spectrum will have a break. The break po- 
sition £br is given by equation ( |82| ) for the idealized case 
of a burst with a constant flux; it appears at ~ 10 MeV 
in the beginning of the GRB and then slowly shifts with 
time to higher energies. In highly variable bursts, future 
time-resolved spectroscopy should show a positive corre- 
lation between ebr and the instantaneous flux. Once the 
break is observed one can evaluate the wind density. 

The main observational effect of the radiation front is the 
strong softening of the early blast wave emission (which would 
otherwise be in the hard X-ray band). Owing to this softening 
the blast radiates in a different spectral window compared to 
the prompt GRB and it can be studied separately in simultane- 
ous observations. Observations in optical - soft X-ray bands at 
early times (less than ~ 1 min) can help to establish the nature 
of the GRB progenitor — as we discussed here a wind from a 
massive progenitor should have clear signatures. 

Early optical emission has already been detected in 
GRB 990123 (Akerlof et al. 1999) and it likely comes from 
the blast wave rather than internal dissipation in the ejecta since 
there is no correlation between the optical light curve and the 
prompt GRB. The optical emission can be produced by the re- 
verse shock in the ejecta (e.g. Sari & Piran 1999). The results of 
the present paper suggest an alternative interpretation: the soft 
emission is produced by the forward shock at its early stage 
when it propagates in the preaccelerated and pair-loaded envi- 
ronment. 

The emission from the blast hardens fast when the blast wave 
propagates in the e ± loaded zone R < R\ oa d and the observer 
will see the whole spectrum from optical to X-ray bands. This 
broad-band emission can overlap with the prompt GRB and its 
X-ray component can affect the measured GRB spectrum. In 
particular, the early external shock may generate the soft X-ray 
excesses detected in GRBs. 



The simple blast wave model constructed in this paper as- 

1/2 — 7/3 

sumes a burst duration t b < t* b = AE^ T ej2 ' s (cf. § 7.2). 
After a redshift correction (1 + z) ~ 3, t* b is still smaller than 
r-j 10 s, while a number of observed GRBs are longer than 10 s. 
The extension to % > t b should be simple since the found front 
structure in terms of £ (§§ 4-6) will hold regardless tj,. How- 
ever, the blast wave will be affected: (1) i? gap will decrease 



(eq.|90|), (2) f A (-R) [and hence j(R)] will change, and (3) the 
sweeping-shell approximation (§ 8) will not work; instead one 
has to deal with the full hydrodynamical problem with the for- 
ward and reverse shocks. 
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APPENDIX 

SCATTERING AND PHOTON-PHOTON ABSORPTION 



Compton scattering 

5 All the observed times given in this section are for GRB redshift z = 0; they should be multiplied by (1 + z) if z > 0. 
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The differential cross-section for Compton scattering (defined in the electron rest frame) is given by (Jauch & Rohrlich 1976) 

V * = ^ + ^-2(l- A1 ) + (l-M) 2 , %= , 1 ~ , (1) 

C Ccr- 



da 



dp 8 CTT 



i - c(i - A) 



where e and e sc are the photon energies before and after scattering respectively (as measured in the electron rest frame), and p is the 
cosine of the scattering angle in the rest frame. In our problem the scattering medium is cold and has a bulk velocity (3 parallel to the 
direction of the primary collimated photons. The rest-frame magnitudes are then related to the lab ones by 



i = j(l-P)e, e sc = 7(1 - Pfj,)e s 



/i = 



1-/V 



da 

dfi dfj, dp, 



ii 8 1 V e ) 1-8 ' 



The Lorentz-invariant total cross-section is given by 



3 <tt 



2_ 2 

e~ e 2 



ln(l + 2e) 



1 



(2) 



(3) 



e 2(l + 2f) 2 _ 
Saturation of radiative acceleration 

When the medium accelerates, the typical photon energy in the medium rest frame is redshifted well below m e c 2 and the scattering 
occurs with Thomson cross-section, da/dp = (3/8)<tt(1 + p 2 ) and okn = ct- With increasing 7 the finite collimation angle of 
the radiation intensity 1(9) becomes important and the efficiency of radiative acceleration drops. The radiative force accelerating the 
electron is (Gurevich & Rumyantsev 1965) 



dp 
dt 



= v 7 2 [Mi + P 2 ) - ih + , h = J 1(6) co & k edn. 



(4) 



Here 9 is the angle between the ray and the radial direction. Note that the net flux F = ii . If the radiation field is perfectly collimated 
(Id = I\ = I2 = F) then dp/dt — (or l c )F(l — (3)/(l + /3). For a finite collimation there exists a frame with a velocity /3 sat where 
the radiation flux vanishes. Assume that radiation is isotropic in this frame and has moments I\ = and I 2 = Iq /3, and compute Ik 
in the lab frame. The Ik are components of the stress-energy tensor of radiation, I — cT 00 , I\ = cT 0x , I 2 = cT xx (the x-axis is 
chosen along the radial direction). From the transformation law, T lk — T lm K\ Aj^ where A is the Lorentz matrix, one gets 



A 



7 2 at 7 0, 



g/3sat7 sa t-' r 0, 



Using these relations and substituting F = Ji we find 



dp 
dt 



(Tt-F 1 7 2 



■(/3sat-/?)(l-/3/?sat) 



1 + /3 



1 



(1 - Aat) 2 



l + (3 



where the approximate equalities make use of 7 sat » 1. The acceleration vanishes when 7 reaches 7 sat . 

7 — 7 absorption 
The 7 — 7 opacity seen by a scattered photon (p,, e sc ) is given by 

K^-y(/i, 6g C ) 



Osat 



(5) 



(6) 



ethr m « C 6 



£thr 



(1 - ^)e sc ' 



where cr 77 is the cross section for 7 — 7 pair production (Jauch & Rohrlich 1976) 



( \ - 3aT 

T ll\ e c) — „ 2 



2 + |-l)m( £c + v ^-T)-(i + 4 




(7) 



(8) 



and e c = (e/ethr) 1 ^ 2 is the energy of the interacting photons in their center-of-momentum frame. The mean energy of the photon 
absorbed by our photon (fi, e sc ) is 



e a bs(M 



e(e c )P(e c )de c = x^thr, P(t c )<ie c = 



F e cr 77 (e c ) dlne c 



f/ ebr/ethr F e a 77 (e c )dlne c 



(9) 



Here P(e c ) is the probability of 7 — 7 interaction with given e c and e(e c ) = ethi-e 2 .- Thus defined numerical factor \ depends on the 
spectrum shape F e . If the absorbing radiation has a power-law spectrum F f = i*\e~ Q one gets at e t hr -C ebr 



_ <i>(a) /e t hr\~ a 
Ax I 2 J ' 



i' OO 

= 2 i- a « 
Ji 



77 -2a-l d 



X = 



<Ka-l) 
20(a) : 



(10) 



20 



RADIATION FRONT OF GAMMA-RAY BURSTS 



where Ai = m e c 3 / Fictt- The numerical factor (j)(a) is with high accuracy (< 0.3% for < a < 6) approximated as <j>{a) 
(7/12)2-"(l + a)" 5 / 3 [Svensson 1987, eq. B6 where 77 = 2 a+1 (a + 2)- 1 0]. It gives \ = (1 + a _1 ) 5/3 . 
The mean energy and momentum of the pair created when the scattered photon gets absorbed are 



e±(fi, e sc ) = (e sc + e abs ) m e c 2 , p±(p, e sc ) = (m £ sc + e a bs) m e c. 



(ID 



In the rest frame of the medium, the mean Lorentz factor and momentum per injected particle are given by Lorentz transformation 
of the energy-momentum vector, 



2TO e c 2 7i n jO, e sc ) = 7 (e± - (3cp±) , 2p inj (/i, e sc ) = 7 (p± - /3^±J . 



(12) 



It is straightforward to further average 7i n j and p; n j over the whole primary spectrum and scattering angles. 

The effective Klein-Nishina cutoff 

To the first approximation, 6kn ~ 7- This estimate is sufficient if ekn 3> 1, far from the spectrum peak. However at 7 ~ 1 we 
have £kn near the peak and the results of the analytical model in § 5 are sensitive to the exact position of £kn- The effective ckn 
depends on what we calculate. In calculations of h+ (eq. || of the paper), da enters in combination with opacity k 77 seen by the 
scattered photon. Let us assume j3 = and compute the average 



f7KNK 77 (e) = / d(TK 77 (/i, e sc ) 



CTT 



(a 2 )e Q2 X Q2 (e), 



X a = 



(i-m)- 



d/i 

1 



d/i 



(13) 



Here bar denotes the averaging over scattering angles. At e < 1 we are in the Thomson regime with X a = X^ — (3/8) J(l — 
/i) a (l + /i 2 )d/i, e.g. Xj — 7/5. The Klein-Nishina correction is important (~ 1/2) already at e ~ 0.1. It is due to two effects: (1) 
the scattering angle is reduced [and k 77 oc (1 — fi) a2 ] and (2) the total cross-section <tkn is reduced. Photons of energy (e, e + de) 
contribute to h+ with approximate weight oc (F e /e)o'KNK^ 7 ", therefore we define the effective ckn for e ± loading as 



F 
e 



Xa 2 — 



e 



'a 2 ' 



£KN = 



(0!2 - a\) 



1 x 

"2 



" 1/(Q2-Ctl) 



0.4, 



(14) 



where f e = F t jF\. We get £kn ~ 0.4 for ax = and ct2 = 1.5, 2, 2.5. 
When calculating P± (eq. n8 of the paper), we need to evaluate 



cr K NK 77 p±(e) = / dcTK 77 (/i,e sc )p ± 



<7t 

a7' 



(a 2 -l)e a2 - 1 X a2 _ 1 (e), 



(15) 



where we neglected the K 77 jie sc term in K 77 p±. Photons of energy (e, e + de) contribute to P± with approximate weight 
oc (F e I e)(7KN k 7 7P± (see eq. jlaof the paper), and the effective exN is given by 



e 



de l£ £ a 2 -l 



^a 2 — 



de— e" 2 " 1 
e 



(a 2 -ax -I) de/ £ e Q 



«2 — 1 



A a 2 -1 



0.4. 



(16) 



We assumed here that a 2 — ax > 1. Again we get ckn ~ 0.4 for ai = and a 2 — 1.5, 2, 2.5. 
In a similar way, one can evaluate the effective ckn for P sc (see eq. [l7| of the paper), 



0kn(1 - A*—) = <?tZ, 



deF f Z 



deF f Z T , 



de/ e 



Z 



0.7. 



(17) 



Here Z T = (1 + /3)- 1 = 1 at /3 = 0. 

When the medium accelerates, ekn increases well above unity and ends up outside the spectrum peak; then its exact value is 
unimportant. A sufficient approximation at /3 > is £kn(/3) = £kn(0)7(1 + (3) i.e. just the Doppler shifted value found at (3 = 0. 



REFERENCES 



Akerlof, C. et al. 1999, Nature, 398, 400 
Blandford, R. D., & McKee, C. F. 1977, MNRAS, 180, 343 
Chevalier, R. A., & Li, Z.-Y. 1999, ApJ, 520, L29 
Dermer, C. D., & Bottcher, M. 2000, ApJ, 534, L155 

Gurevich, L. E., & Rumyantsev, A. A., 1965, Sov. Physics - JETP, 20, 1233 
Jauch, J. M., & Rohrlich, F. 1976, The Theory of Photons and Electrons, 

Springer, New York 
Madau. P., & Thompson, C. 2000, ApJ, 534, 239 
Meszaros, P., Ramirez-Ruiz, E., & Rees, M. J. 2001, ApJ, 554, 660 



Piran, T. 1999, Phys. Rep., 314, 575 

Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., Paciesas, W. S.. 

& Band, D. L., 2000, ApJS, 126. 19 
Rees, M. J., & Meszaros, P. 1992, MNRAS, 258, 41 
Sari, R., & Piran, T. 1999, ApJ, 517, L109 
Smolsky, M. V, & Usov, V. V. 2000, 531, 764 
Svensson, R. 1987, MNRAS, 227, 403 
Thompson, C, & Madau, P. 2000, ApJ, 538, 105 (TM) 
Woosley, S. E. 1993, ApJ, 405, 273 



